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Abstract 

The two-point correlation function in tliin films is studied near the critical 
point of the corresponding bulk system. Based on fieldtheoretic renormaliza- 
tion group theory the dependences of this correlation function on the lateral 
momentum, the two distances normal to the free surfaces, temperature, and 
film thickness are determined. The corresponding scattering cross section 
of X-rays and neutrons under grazing incidence is calculated. It reveals the 
various singularities of the two-point correlation function. 
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I. INTRODUCTION 



Structural properties of condensed matter depend sensitively on the space dimension d. 
Thin films offer the opportunity to reveal this dependence. By varying the film thickness 
L one can interpolate smoothly between d = 2 and d = 3. For crystalline materials this 
variation can be accomplished with atomic resolution by using molecular beam epitaxy [|T|. 
As an alternative, which is applicable also for fluids, thin films can be built up via wetting 
phenomena where the film thickness is controlled by temperature or chemical potentials . 
Once such films are prepared the dependence of their structural properties on the space 
dimension can be studied particularly clearly close to phase transitions. For first-order 
phase transitions the main influence of a variation of the film thickness is to shift the phase 
boundaries in the phase diagram (see, e.g., capillary condensation ||^ or the shift of the 
melting curve [Q) without changing much the /oca/ properties of condensed matter. In rare 
cases, however, even the character of the phase transition can change as function of L; see, 
e.g., the possibility of continuous melting in = 2 as opposed to (i = 3 or the crossover 
from a first-order phase transition in = 3 to a second-order phase transition in c? = 2 at a 
certain thickness of a slab of the 3-states Potts model 0. 

In the case of first- order phase transitions the robustness of the local structural proper- 
ties with respect to changes of the film thickness is due to the smallness of the correlation 
lengths which characterize these systems and — putting aside possible wetting phenomena 
— thus severely limit the propagation of the structural changes, which necessarily occur near 
the confining surfaces of the film, into the interior of the films. In contrast, second- order 
phase transitions are characterized by diverging correlation lengths which affect not only 
the location of phase boundaries but in addition lead to pronounced changes in the local 
properties even deep in the interior of the films if the critical point is approached. These 
effects are thus not only particularly suitable to shed light on the aforementioned depen- 
dence of the structural properties on space dimension but they offer an additional advantage: 
the divergence of the correlation length as function of temperature upon approaching the 
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critical point leads to universal behavior which makes a quantitative comparison between 
theoretical predictions and experiments much easier as compared with systems exhibiting 
first-order phase transitions which are characterized be several competing length scales of 
comparable, atomic size which are difficult to determine accurately and to vary systemati- 
cally and independently. 

A sizable body of theoretical research has emerged describing continuous phase transi- 
tions in thin films (see, e.g., Refs. HtHT^)- Initiated by the theory of finite-size scaling (see, 
e.g., Refs. [|T^ -p3|), inter alia the shift Tc{L) of the critical temperature with respect to its 
bulk value Tc = Tc{L = oo) [^^-|26|, the magnetization [^,^ as well as the free energy. 



the Casimir force, and the specific heat |29-35| have been analyzed. Here we emphasize 



that in order to observe universal film behavior the thicknesses L of the films have still to 
be large on an atomic scale. This is assumed to be the case throughout our analysis. The 
analytic description of the dimensional crossover between d = 3 critical behavior near Tc 
and the d = 2 critical behavior near Tc{L) poses still a challenge [0,0 which has not yet 
been overcome with satisfactory quantitative accuracy. Numerous experiments (see, e.g., 
Refs. PB|-^) and simulations (see, e.g. Refs. P^-^) have been carried out to test these 



theoretical predictions. They lend support to the finite size scaling theory but still pose a 
puzzle as far as detailed quantitative agreement is concerned. 

The vast majority of these studies is devoted to integral or excess quantities without spa- 
tial resolution. However, the studies of /oca/ critical properties, such as of one- and two-point 
correlation functions, near a single surface have revealed a wealth of universal phenomena 
featuring numerous surface critical exponents and interesting crossover phenomena — on the 



scale of the bulk correlation length ^ — between surface and bulk critical behavior p7| , ^ 
the integral and excess quantities offer either no or only very limited access to these local 
properties. 

The successful development of surface specific X-ray and neutron scattering techniques 
based on exploiting total external reflection at grazing incidence has proven to be very 
fruitful, inter alia, for facilitating the quantitative comparison between experiments and 



theoretical predictions of the local critical behavior near interfaces ||49| , |50[] . These scattering 
techniques allow one to determine order parameter profiles normal to the surface and the 
depth resolved lateral two-point correlation function. In the present context such experi- 



ments have been carried out successfully for the binary alloy Fe^Al |]5T|-p^ and, by using 
truncation rod scattering, for FeCo which exhibit continuous order- disorder transitions 
in the bulk. In the case of Fe^Al the cusplike surface singularities of the momentum and 
temperature dependence of the two-point correlation function turned out to be in excellent 
agreement with the theoretical predictions ||55|j56[| . The fact, that due to the occurrence of 
surface segregation suitable choices for the crystallographic orientation of the surface allows 
one to switch between the different surface universality classes corresponding to free bound- 
ary conditions and boundary conditions with surface fields, respectively, of the same bulk 
sample ||57|j58[| , offers wide ranges of interesting comparative studies. 



In view of these developments and in view of the increasing availability of powerful 
synchrotron and neutron sources it appears promising to extend these studies of local critical 
properties to thin films. There are several predictions concerning the behavior of one-point 
correlation functions in thin films such as order-parameter profiles P7| , |28|P5|] and energy- 



density profiles However, on the level of the two-point correlation function so far only 
very little is known. This function depends on the lateral distance xy = Xy^' — x.\^^^ between 
the two points xi = {x.\^^\ Zi) and X2 = (X||^'',2;2) (or equivalently the lateral momentum 
p corresponding to the d — 1 translationally invariant directions), the coordinates Zi and 
Z2 perpendicular to the parallel surfaces of the film, the film thickness L, and temperature 
t = (T — Tc)/Tc (or equivalently the bulk correlation length ^ = ^o^~'^)- Since a full sweep 
of this large parameter space is practically not possible for computer simulations, we have 
applied fieldtheoretic techniques which provide analytic access to the full parameter space. 
This approach encompasses nonperturbative features such as scaling properties and short- 
distance expansions as well as an explicit and systematic perturbative result to first order 
in e = 4 — d. The latter serves as to corroborate the nonperturbative results and to provide 
numerical results which are not accessible by general arguments. These explicit calculations 



are carried out for the fixed point of the so-called ordinary transition for both confining 
surfaces in the classification scheme of surface critical phenomena corresponding to free 
boundary conditions on both sides. This is applicable to thin antiferromagnetic films near 
their Neel temperature, to ferromagnetic films near their Curie temperature in the absence of 
external bulk and surface fields, and to thin films of binary alloys near their continuous order- 
disorder transitions. Among the numerous order-disorder phase transitions in binary alloys 
only a few are of second order including Fe^Al |]60H6^, FeCo [^, CuZn [^, and FeAl . 



Both the B2 — DO3 transition in Fe^Al and the A2 — B2 transitions in FeCo, CuZn, and 
FeAl belong to the Ising universality class . For the A2 — B2 transitions it is predicted 
theoretically that the (110) surface belongs to the surface universality class of the ordinary 
transition whereas the (100) surface exhibits the so-called normal transition associated with 



the presence of an effective surface field [57,58|. Indeed truncation rod scattering at the 



FeCo (100) surface has provided clear evidence for the presence of an effective surface field 



5J] above T^. although the expected associated crossover from ordinary to normal critical 
behavior could not yet been resolved experimentally in an unequivocal way. The results 
of the diffuse scattering of X-rays under grazing incidence from the (llO) surface (equivalent 



to the (110) surface) of Fe^Al ||^ are in excellent agreement with the theoretical predictions 



55| , |56[| for the ordinary transition. But even for Fe^Al (110) a residual order parameter 



above has been reported [pl| , p2| . Thus it still remains to be seen theoretically whether 
for the B2 — DO3 transition in Fe^Al, in contrast to the A2 — B2 transition in FeCo, the 
(llO) surface can support a weak effective surface field. In view of this state of affairs our 
present result are expected to be closely applicable to thin films of Fe^Al, FeCo, CuZn, and 
FeAl bounded by (110) surfaces on both sides. Among them Fe^Al and FeCo appear to be 
the most promising candidates because the others exhibit strong surface segregation. For 
an assessment of the possibilities to probe critical magnetic surface transitions by grazing 



incidence of neutrons see Ref. 



In view of the aforementioned difficulties concerning the analytic description of the di- 
mensional crossover we confine our analysis to the temperature range T > T^. We note that 



elements of the perturbation theory for thin critical films can be found in Ref . |p8| . We had, 



however, to carry out our own approach because the representation given in Ref. |68| is not 



suited for making predictions for the scattering experiments and because Ref. contains 
errors. Finally we note that experience tells that calculations carried out for the spherical 
model, as have been done for the present system , lack the quantitative reliability needed 
for comparison with experiments and simulations. 

In order to encourage future scattering experiments for critical thin films and to facilitate 
an explicit quantitative comparison of such data with the present theoretical predictions, 
we have calculated the singular contributions to the scattering cross section for X-ray and 
neutron scattering under the condition of grazing incidence based on our results for the 
critical two-point correlation function in thin films. This allows us to describe the conditions 
under which the various singularities of the two-point correlation function become visible in 
scattering data. 

This introduction is followed be three sections, the Summary, and four Appendices. 
In Sec. we introduce the fieldtheoretical model. The two-point correlation function is 
discussed in Sec. |TT| and in Sec. we investigate the scattering cross section. Relations 
between bulk and film amplitudes are derived in Appendix 0, explicit one-loop results 
are presented in Appendix and Appendices and contain details required for the 
calculation of the scattering cross section. 

II. FIELD-THEORETICAL MODEL 

The leading critical behavior in a film follows from the statistical weight exp(— 7^{$}) for 
the configuration $(x) = (0a(x), a = 1, . . . , n) of a n-component field, which is proportional 



to the order parameter, where p8|j30| , |32[] 



7t:{$(x)} = J d'^-^xii dz Q(V$)2 + + |($2)2 _ h ■ $j (2.1) 
+ j d'^-^x\\ [^"^^z = 0) - hi ■ $(2 = 0) + ^<l>\z = L) - hi ■ $(2 = L)) 
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with space dimension d and position vector x = (xy, z) of d—1 parallel and one perpendicular 
components. The z integration extends over the interval [0, L], where z = and z = L give 
the positions of the film surfaces, r is the temperature parameter such that in the bulk r = 
marks the transition temperature within mean-field theory. The coupling constant g > 
ensures the stability of the statistical weight below the transition temperature, i.e., for r < 0. 
c denotes the surface enhancement, h and hi are bulk and surface fields, respectively. We 
focus on the ordinary transition at zero fields, i.e., we adopt the fixed point value c = oo for 
the surface enhancement and set h = hi = 0. After carrying out a Fourier transformation 
with respect to the d — 1 directions exhibiting translational invariance parallel to the surfaces 
the mean-field propagator for the disordered phase (r > 0) in p-2;-representation is given by 
4|,|70l 

Gn{p,z,,Z2,L,T)= I d^-^x. e^P-^ii (<l>(x,|, ;2i)<l>(0, Za)) (2.2) 



}_l -b\zi-Z2\ _ p-b{zi+Z2) 

2b 



+ T^KT—. ' ^ = VV+^. 



Q~b(zi-Z2) _j_ g-b(22-2l) _ Q-b(z-i+Z2) _ Qb(zi+Z2) 

e^^L - 1 

The first exponential function corresponds to the bulk part followed by the contribution from 
the surface at z = 0. Both exponentials together give the propagator for the ordinary transi- 
tion of the semi-infinite system (L = oo). The remaining ratio carries the L dependence. The 
propagator satisfies the Dirichlet boundary conditions Gd{z = 0) = = G^iz = L). Equa- 
tion ( p.2|) represents the mean-field approximation for the two-point correlation function in 
the film corresponding to the critical behavior in ci = 4. The non-Gaussian fluctuations in 
d = 3 are taken into account approximately by the one-loop correction which amounts to 
the first term in a systematic expansion in terms oi e = 4 — d: 

GbareiP, Zl, Z2, L, T, q) = Gd{P, Zi, Z2, L, t) (2.3) 



gn + 2 f d'^-^q 



L 



2 3 J (27r)^-i J GDip,Zi,z,L,T) Goiq, z, z, L,t) Goip, z, Z2, L,t) + 0{g^). 


As regularization scheme we use dimensional regularization by analytic continuation in the 
space dimension d = 4 — e. As long as Zi and Z2 are both off the surfaces only bulk 



singularities occur. We absorb the corresponding poles in e by minimal subtraction through 
the standard Z factors: 

= V"", 9 = /i^2'^7r'^/2Z„M, r = fx^Z^t, (2.4) 



where n is the momentum scale and the bulk Z-factors are |7T| 



Z^ = l + 0{u'), Z^ = l + ^^- + 0{u'), Z, = l + ^^ + 0{u'). (2.5) 

O 6 o 6 



The renormalized correlation function reads (see Eq. ( |2.4| )) 

G{p, Zi, Z2, L, t, U; /i) = Z^^GbareiP, ^2, L, T, g) (2.6) 

which is valid in all orders of perturbation theory. The solution of the corresponding renor- 
malization group equation leads to the following scaling property: 

G{p,z,,Z2,L,t;fi) = g^p-'+'^gM,zi/^,Z2/^,L/0. (2.7) 

This holds at the fixed point u* = + C(e^) and involves the bulk correlation length 
^ = ^Qt"'^, the exponents r] = C(e^), and z/ = \ + \^^^+0{e'^). With suitable normalization 
(see, c.f., Eq. (p.l3|) ) the scaling function gi is universal. The amplitude Qi, which is fixed by 
this normalization, and the amplitude carry the nonuniversal scaling factors. We fix by 
defining ^ as the so-called true correlation length |]72[ so that = yU~^(l + i^7i]^|(l ~ Ce)^ + 
0{e^)). This expression for allows one to express the momentum scale /i introduced in 
Eq. (p.4|) in terms of the experimentally accessible, nonuniversal amplitude 

= (eo^)-^(l + - CE)e + 0{e^)). (2.8) 

If subsequent formulae contain the momentum scale /z explicitly it is to be replaced by Eq. 
( p.8|) ; moreover we omit /i from the explicit list of variables of G. 

Depending on the problem under consideration it is often advantageous to use different 
but equivalent representations of the correlation function such as 

G{p,zi,Z2,L,t) = guzl''^guipz2,zi/^,Z2/^,Z2/L), (2.9) 



Gip, zi, Z2, L, t) = guiL^-'^guiipzi, zi/L, Z2/L, (2.10) 

G{p,Zi,Z2,L,t) = gjy^^-''gjy{pL,pzi,pz2,^/L), (2.11) 

and 

G{p,z,,Z2,L,t) = g^p-^+^g^{pi,p{zi - Z2),p{zi + ^2), VO- (2-12) 

The nonuniversal amplitudes and the universal scaling functions Qx-, x = I , II , III , IV, V, 
are fixed by the following normalizations: 

lim lim \im gi{a, j3,'j = P,S) = 1, (2-13) 

a— >oo /3— >oo 5— »oo 



lim lim lim gu{a, /?, 7 = 5) =: (7„(0, 0, 0, 0) = 1, (2.14) 

a— »0 /3^0 5^0 



lim lim (?:„(«, /3= 1/2,7 = /3 = 1/2,5) = 1, (2-15) 



lim lim gj^{a, /3 = a/2, 7 = ^ = a/2, 5) = 1, (2.16) 

and 



lim lim lim lim gy{a, ^,^,6) =: gv{oo,0, 00,00) = 1. (2-17) 
The universal scaling functions g^ can be expressed in terms of each other because in Eqs. 



and (|2.9| ) - (|2.12|) the left hand side is the same quantity and the sets of scaling variables 
are complete, i.e., from each set one can form any of the others by a suitable combination 
of variables. 

Since the nonuniversal amplitudes Qx correspond to the same correlation function 
G{p, Z\, Z2, L, t) and because the scaling functions fixed by the normalizations in Eqs. ( p.l3|) 
- ( |2.17| ) are universal, their ratios Qx/Qx' are universal numbers. Thus the knowledge of one 
of them and of the corresponding universal scaling functions determines all the others. 
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Moreover, as discussed in Appendix 0, all nonuniversal amplitudes Qx are determined 
by any pair of nonuniversal scale factors which characterize the critical bulk properties. A 
transparent and experimentally directly accessible choice for the latter is the nonuniversal 
amplitude B of the leading temperature singularity of the field (0(x)) in the bulk below Tc, 

(</>(x)) = B{-tr, (2.18) 

and the amplitude of the true correlation length above Tc. In terms of these quantities 
one has 

= B^^+y-'^+^U (2.19) 

where W is a universal number, whose value U ^ 1.58 is derived in Appendix ^ based on 
Eq. ( |2.17| ). In the following most of our analysis focuses on the scaling function used in 
Eq. (|2.9|) . For that case one finds (see Appendix ^ the universal ratio 

g,,/g, = 2(1 + + 0{e')). (2.20) 

With these results we finally obtain 

G{p,z,,Z2,L,t) = B\itY~'n{z,/itY-^g,,{pz2,z,/^,Z2/i.Z2/L) (2.21) 

where TZ = 2u(^l + + 0{e^)^ ^ 4.21 is a universal number. Thus in all our subse- 
quent formulae for film properties their absolute values are determined and fixed by the two 
nonuniversal bulk amplitudes B and ^q. 

The actual order parameter OP for a particular second order phase transition is propor- 
tional to the field (p introduced in Eq. ( p.l| ), i.e., QP(x) = &0(x). The value of b depends 
on the particular system (binary alloy, liquid, ferromagnet etc.). Moreover, any rescaling 
of 6 by a dimensionless number renders another order parameter OP which is equally valid 
for describing the singular behavior of the phase transition. We emphasize that Eqs. ( p.9[ ), 
( Prr^ ), and ( pUD remain vahd if G is replaced by (OP(x)QP(x')), (0(x)) by (OP(x)), and 



B hj B' = bB; these replacements have to be carried out if the present fieldtheoretic results 
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are used to interpret, e.g., the intensity of scattered X-rays or neutrons (see, c.f.. Sec. |rv| ). 
The actual choice of the QP, as it enters into the expression for the scattering cross section, 
is borne out and tight to the relation (QP(x)) = t)^. 



III. EXPLICIT PROPERTIES OF THE TWO-POINT CORRELATION 

FUNCTION 

The discussion of the correlation function consists of three parts. First, we set zi = 
and analyze its nonanalytic behavior in certain limits. Then, we take into account the 
case z\ 7^ 2:2, which serves to understand the correlations perpendicular to the surfaces. 
Moreover, the discussion of this latter case turns out to be very useful for carrying out the 



integrations appearing in the scattering cross section to be analyzed in Sec. The film 
excess susceptibility is discussed in the last part. 



A. Lateral two-point correlation function for z\ = Z2 

In order to investigate various asymptotic properties of the lateral behavior of the two- 
point correlation function we resort to short distance expansions (SDE) |j73|, distant wall 
corrections (DWC) and results of the perturbation theory supported by appropriate 
exponentiations of the explicit e-expansion results. With z\ = Z2 = z, in the present context 
a representation of the form 

G{p, z, L, t) = gnz'-''>g{pz, zli, zjL) (3.1) 

is useful. According to Eq. ( |2.9| ) one has g{u,v,w) = gu{u,v,v,w) with (7(0, 0,0) = 1 (Eq. 
( p.l4| )). For semi-infinite systems, i.e., L = 00 the SDE in the cases t = 0, p — > and p = 0, 



t ^ [[M| , [r5[ leads to the asymptotic behaviors 



G{p, z,L = oo,t = 0) = Gnz^'^giiu = pz) (3.2) 

C„^i-^[l + Ai(p^)-i+''ii + ...] 
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and 



G{p = 0,z,L = oo,t) = G,,z^-^g2{v = z/i) (3.3) 

= Guz'-^l + B,{z/^+)-^+^H-'-- + ...], 

respectively, with 711 = i/(?7|| - 1), gi{u) = g{u,v = 0,w = 0), gi{0) = 1, g2{v) = g{u = 
0,v,w = 0), and (72(0) = 1. In the case p = 0, t = one has 

G{p = 0, L, t = 0) = Guz'-^'g-siw = z/L) (3.4) 

with gz{w) = g{u = 0,v = 0,w) so that gs{0) = 1. In order to infer the first nontrivial 



dependence on L for L ^ 00, according to Eq. (|3.4| ) one can equally consider the limit 
z ^ for L fixed. To this end we consider the SDE of the renormalized film correlation 
function in real space: 

(</)(x||,z)0(O,z))-/i-2(/x^)2(--^')(0^(x|,)0^(O)) (3.5) 
= /'2(/i2)'("=-")(/iX||)-2^'=r(x||/L). 

Here 0_l denotes the normal derivate of taken at one of the surfaces and Y{y) is a dimen- 
sionless scaling function for the film which is universal up to nonuniversal prefactor. The 
scahng dimensions of (p and (f)± are x = ^{d — 2 + rj) and Xg = ^{d — 2 + m) , respectively. 
The scaling function Y{x\\/L) describes the influence of the distant wall at z = L on the 
lateral correlations close to the near wall at 2; = 0. In order to obtain its leading asymptotic 
behavior for xy/L — >• we use the identity 

G{p = 0,z,L,t = 0) = G{p = 0,z,L = oo,t = 0) (3.6) 

The first term on the rhs is equal to GuZ^~^ (compare Eqs. ( p.2| ) and (pl3|)). The leading 
correction is given by using the SDE in Eq. (|3.5|) for the second term: 
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dG{p = 0,z,L',t = 0) ^^, 
dL' 



d 



(3.7) 



0) = 1 + CiW-^+i\\ where 



with C = d'^-^y y-^'^---+m)Y'{y). Thus we find g^{w 

II is a universal number, i.e., 

G{p = Q,z,L^oo,t = Q)=g,,z^-'^[l+Ci{z/L)-^+'^\\ + ..]. 



(3.8) 



Finally we note that due to the normalization (/(O, 0, 0) = 1 the scaling function g{u, v, w) is 
given by the ratio g{pz, z/C^, z/L) = G{p, z, L,t)/G{p = 0, z, L = oo,t = 0) from which the 
prefactors QuZ^~^ appearing in Eq. (|3.1| ) drop out. The e-expansions of the amplitudes of 
the leading asymptotic terms follow from Eqs. ( |B6|) , (p7|), and (p8|)) in Appendix p 1| : 

2 







1 + e 


Bi = 




1 + e 


Ci = 




1 + e 



n 



'l-CE-\n2) + 0{e' 



n + 

n + 2 /tt^ 



(3.9) 



Ce + 2(^2 + /i) - 1 + 0{t 



n + 8 V 18 



0.5772 is Euler's constant, S2 ^ 0.083 and h ~ 0.287 are given by Eq. in 



Appendix P 1| . Within the e-expansion the full forms of the scaling functions gi{u), g2{v), 
and gsi^w) can be found in Appendix P 1| (see Eqs. ( pi]) - (p3D). 

In Fig. |l] we display the three scaling functions gi, i = 1,2,3, (Eqs. ( pT]) - (p3|)) 



corresponding to Eqs. (p.2|), ( p.3|) , and ( |3.4| ) as obtained within mean- field theory (MET), 
i.e., for e = and from renormalization group guided perturbation theory (PT) as well 
as their leading behavior gi{xi ^ 0) = gi,i{xi), Xi = u,X2 = f , ^3 = w. Within MET 
the three scaling functions have the same limiting form for small scaling variables with 
Ai = Bi = Gi = —1 and the critical exponent t]\\ = 2. Beyond MET, in Eig. |^ we 
use ?7|| = 1.48 as the best available estimate whereas the amplitudes are evaluated in 

13 



first order in e (Eq. for (n, e) = (1,1)) so that Ai ^ -0.9099, Bi ~ -1.1409, and 

Ci ~ —0.9035. Within mean-field theory gi = g2 and the leading asymptotic behavior g-^ i 
provides already the full scaling function g^. Beyond MFT there is a small difference between 
(73 and g^^. This difference is much bigger for the scaling functions gi and g2 describing the 
semi-infinite system. 

The above discussion demonstrates that, for z fixed, the two-point correlation function 
G{p, z, L,t) has a finite value G{p = 0,z,L = oo,t = 0) which is attained via cusplike 
singularities: ~ j9"^+''ii {p 0, 1/L = 0,t = 0), ~ (l/Ly^+^w (1/L ^ 0,p = 0,t = 0), 
~ (l/,^)~^'^''ii {t — > 0,p = 0, 1/L = 0). In terms of these variables the critical exponent is 
the same for all three cases and only the amplitudes differ. These singularities remain if 
only one out of the above three variables is zero and the remaining two both vanish. This 
behavior, which includes the smooth interpolation between the corresponding amplitudes, 
is described by the scaling functions hi, h2, and of two variables instead of the scaling 
functions with one variable as gi{u), (72 (f), and (73 (w): 

G{p, z, L = 00, t) = QijZ^^'^hi{u,v), hi{u,v) = g{u,v,w = 0), (3.10) 

G{p, z, L,t = 0) = QuZ^~'^h2{u,w), h2{u,w) = g{u,v = 0,w), (3-11) 

and 

G{p = 0, z, L,t) = QjiZ^~^h3{v,w), h^i^VjW) = g{u = 0,v,w) (3.12) 

with u = pz, V = z/C,, and w = z/L. All three scaling function can be obtained from Eq. 
( |B17| ). Since the discussion of all three scaling functions is analogous we demonstrate our 
analysis only for h^iy^w). We introduce polar coordinates uj and ip 

UJ = Vv"^ + w'^ = z\J -|- L^2^ = arctan(t>/w) = arctan(L/^) , (3.13) 
v = to" sin yj, w = UJ cos 

which leads to 
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hsiv^w) = hsiuj sin ip,uj cos ifi) = hf]i^^{uj , ip) . (3.14) 

Since the limit u ^ 0, i.e., 1/.^ ^ and 1/L 0, is equivalent to the limit z ^ for ^ and 
L fixed the resulting singularity is compatible with the SDE so that 

hfLi^ - 0, y.) = i/f (v.) + iff + . . . . (3.15) 

The explicit form of the scaling function h^lvjw) as obtained from perturbation theory in 
(9(e) is in accordance with Eq. ( |3.15| ) and renders explicit results for the coefficients Hq^\(p) 
and H[^\vy. 

^o(^) = h'^Xriu = 0, ^) = h{v = 0,w = 0) = l (3.16) 

is independent of ip and equal to 1 due to the normalization g{u = 0,v = 0,w = 0) = 1. 
With this result the e-expansion of iff ^ {(p) follows by comparing the e-expansion of the rhs 
of Eq. (|3.15| ) with the limit — * of the e-expansion of hpJi^^{iu, cp). As expected one finds 
that h[^\p) interpolates smoothly between the value iff "*(</? = 0) = Ci (see Eq. (|0|)) 
corresponding to the amplitude of the singularity ~ (l/L)^^^''" for u = and v = and the 
value Hi {p = 7r/2) = Bi (see Eq. ( p.9| )) corresponding to the amplitude of the singularity 
~ (l/^)^^^''ii for u = and w = 0. In Fig. ^ all three amplitude functions H[^\p), 
Hf\p), and Hf\p) (see Eqs. ( pToD - (|BT^ )) are shown in mean-field theory (MET) and 



in first order in e (PT). Within MET H^\p) of the semi-infinite system is constant and 
nf^iyp) = h[^\p) exhibit a nontrivial dependence on p. Beyond MET all three functions 
interpolate between the amplitudes Ai, Bi, and Ci (see Eq. (|3l^ )) in a nontrivial way. 

In Eigs. 0, H, and|| we display the full scaling functions hi{u,v), h2{u,w), and h^{y,w), 
respectively. In order to obtain such a scaling function beyond the leading asymptotic 
form we first subtract its leading contribution in its e-expanded form in 0{e) from the full 
expression of the scaling function and add the leading exponentiated contribution afterwards. 
This exponentiation scheme is consistent with the explicit expanded form up to and including 
0{t). In Eigs. ^, 1^, and ^ we show cross sections of the three-dimensional plots in order to 
illustrate the emergence of the p~'^^^\\ cusplike singularity upon varying t or L, the (l/^)~^^''ii 
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cusplike singularity upon varying p or L, and the (l/L) ^^''ii cusplike singularity upon 
varying t or p, respectively. 



B. Perpendicular correlations 

In a semi-infinite system the perpendicular correlations in real space define the exponent 
'7± = (^7 + ^||)/2 through the limit G{x\\,Zi —>■ oo,Z2,L = cxD,t) ~ 2;^ ^j^j^ 
Zi fixed. A Fourier transformation leads to the relation G{j) = 0,Zi,Z2,L = 00, t) ~ zl~^^ 
with Z2 fixed and zi 00. Note that in real space G{x\\, zi, Z2, L = 00, t = 0) increases as 
function of Zi for Z2 and X|| fixed, reaches a maximum at a certain value 2;* = 2;2/(x||/2;2) and 
finally vanishes for zi ^ 00. This increase for < zi < zl leads to the divergence ~ zl^^^ , 
1 — r7_L ~ 0.25, of G{p = 0, zi, Z2, L = oo,t = 0) = f dx\\ (0(0, zi)(f){x\\, ^2))- The coordinates 
Zl and ^2 can be interchanged. Actually conformal invariance fixes completely the functional 
form of G{p = 0, Zi, Z2, L = oo,t = 0) (see Ref. [Q). SDE leads up to a constant amplitude 
to the expression 

Gip = 0,Zi,Z2,L = oo,t = 0)r^{ziZ2y-^(Q{zi-Z2)- + Q{z2-zi)^) ' +... (3.17) 

V Zl Z2J 

(see Eq. (4.68) in Ref. |7^). The explicit calculation to first order in e gives 



111-1 



G{p = ^,zi,Z2,L = oo,t = {)) = Q,lziZ2)^[Q{zi-Z2)-^Q{z2-Zi)-\ " (3.18) 

V Zl Z2' 

for arbitrary Zi and Z2- This perturbation theory guided result for d = 3 has a structure 



similar to the exact result from conformal theory in ci = 2 (see Refs. ||76[ and fT^). Therefore 
one is led to the conclusion that Eq. ( p.l8| ) is a good approximation for the exact correlation 
function in c? = 3. Guided by these considerations we find that in the case that the variables 
p, t, and l/L are small but nonzero the explicit results for G obtained from the e-expansion 
can be cast into the following forms: 



^(p^O, zi,Z2,L = oo,t = 0) (3.19) 

. Zl ■ 



^n(e(^2 - Zl)zr^(^)'"''" (1 + Ai(p^2)-'+'"l + . . . ) 
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■vi\ 



and 



G{p = 0, Zi, Z2, L = oo,t ^ 0) 

^ Z\ . 



Z2 



+ . . . 



+ Q{zi - Z2)Z. 



Z2 



l-r)j_ 



1 +5i 



^1 



+ 



(3.20) 



G(p = 0, zi, Z2, L ^ oo,t = 0) 



(3.21) 



6^„(e(z2-^i)^r"(- 

.2:1 



-1- 



+ . . . 



+ 6(^1-^2)4""^ 
^Z2 



These expressions are valid for arbitrary zi and Z2 as long as the scaling variables pzi,2, 
-21,2/^, and Zi^2/L are small. The explicit e-expansion provides the amplitudes Ai, Bi, and 
Ci given by Eq. (^.91). For the special case Zi = Z2 Eqs. ( |3.19| ) - ( |3.21| ) reduce to Eqs. (^.2p, 
(^, and (|X|). In the limits p = 0, t = 0, and L = 00 Eqs. ( PTTUD - QX^ ) reduce to Eq. 
( CT ) (recall r^^ = 

Finally we note that Eqs. (|3.19|) , (|3.20| ), ( |3.21|) , and the full film correlation function 
G{p, Zi,Z2, L, t) up to first order in e (see Eq. ( P17|) in Appendix |B 21) satisfy the so-called 
product rule derived by Parry and Swain for the correlation function algebra of inhomoge- 



neous fluids (see Eq. (2.20) in Ref. |78|): 



G{p, zi, Z2, L, t)G{p, Z2, Z3, L, t) = G{p, Z2, Z2, L, t)G{p, zi, Z3, L, 



t) 



(3.22) 



for all < Zi < ^2 < 23 < L. The second identity derived by Parry and Swain (see Eq. 



(2.21) in Ref. [|78[) is trivially fulfilled in the disordered phase considered here because it 
involves the derivative of the order-parameter profile which vanishes above T^. A nontrivial 
test of this relation would require results for the ordered phase below T^. 
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C. The susceptibility 



As it has become apparent in the previous subsection, the full dependence of the cor- 
relation function G on all its variables p, zi^ Z2i L, and t is rather complicated. Therefore 
it increases the transparency to consider a spatially averaged quantity which still displays 
interesting specific properties of the critical behavior in a film geometry. The singular part 
of the total susceptibility per area defined as 



rL rL 

x{L,t) = / dzi / dz2 G{p = 0,Zi,Z2,L,t) 
Jo Jo 



(3.23) 



provides such a reduced but still interesting quantity in that it depends only on two variables 
L and t. In addition this susceptibility is directly accessible in an experiment which probes 
the response of a thin magnetic film on the applied external field in the limit of vanishing 
field strength. 

From the scaling properties for G one obtains the following scaling property for x (see 
Eqs. ( g:2lD and (|ATID ): 



(3.24) 



where 



fiy) 



dxi / dx2xl ''gu{0,xiy,X2y,X2) 

^0 



(3.25) 



is a universal scaling function. For y oo, i.e., L ^ oo and t fixed the scaling function 
f{y) vanishes as follows: 



fiy^oo)=Ay-'^'^ + By-'+'^ + Cy 



-3+V 



+ 0{e 



(3.26) 



with 



B = -2(^l+i 
C = 4(l-i 



7l) 



(3.27) 



so that with 7 = z/(2 — rj) and 7s = 7 + z/ 

X{L ^oo,t)= B^i^^f+'nl^At-' + Bt-'^ \l + §e-^/« + 0(e^2^/«)j }. (3.28) 

I. 4q L £5 J J 

The first term (~ t^'^) corresponds to the bulk contribution of the total susceptibility. 
(We recall that x is the total susceptibility per area A\\ of one surface and that the total 
volume of the system is A\\L.) The universal amplitude A (Eq. ( |3.27D ) is in accordance 
with the corresponding known universal amplitude ratios [^,^. The second term (~ t^'^^) 
corresponds to the sum of the excess susceptibilities of two semi-infinite systems within the 
surface universality class of the ordinary transition resembling the two bounding surfaces of 
the film. The corresponding universal amplitude B (Eq. (|3.27|) ) of the semi-infinite systems 
is in accordance with the corresponding result in Ref. ||81[. Finally, the last term ~ e~^^^ 



in Eq. (|3.28|) is the actual finite size contribution induced by the finite distance L between 



the two surfaces confining the film. It is interesting to note that the structure of this finite 
size term Ct"'''" exp(— L/^) differs from its counterparts for the free energy and specific heat 
in two respects (see Eqs. (4.8) and (6.14) in Ref. [30(a)]): (i) For ordinary - ordinary 
boundary conditions the finite size terms of the latter two both vanish ~ exp{—2y) for large 
y = L/^. (ii) The prefactor Ct~"''' is replaced by C't~'^y'^^^'^''^ with k = — 2 (free energy) 
and K = as = a + u (specific heat), respectively. From the explicit result in 0{e) we infer 
that in the case of the excess susceptibility this power law in front of the exponential is 
either missing or has an exponent of O(e^). In order to render the comparison between the 
finite size scaling of the free energy and specific heat on one hand and of the susceptibility 
on the other hand more transparent we rewrite the susceptibility as 

where (2 — 77 = 7/z/, 3 — 77 = 7^/ z/) 



xiL, t) = B'iar^'n[j^y^'^{Ay-^/'' + By-'^'^ + g{y)] (3.29) 



g{y) = f{y) - Ay'^^^^ - By~^^\ (3.30) 

The finite size scaling for the singular part Tsmg of the free energy of a film has the similar 
form (d = (2 - a)/z/, - 1 = (2 - as)lv) (see Eq. (4.11) in Ref. [30(a)]) 
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A II is the area of the cross section of the film. Both in Eq. ( p.29|) and Eq. ( p.31|) the first 



two terms correspond to the bulk and surface contribution, respectively. In both cases the 
curly bracket represents a universal scaling function. For the susceptibility the finite size 
part vanishes as 

g{y ^ oo) = Cy-"''/''e-^ + 0{e-^y) (3.32) 
whereas for the free energy one has 

@{y ^ oo) = c'y-(-=-2)/(2-)e-2i/ + 0{e~^y). (3.33) 

At this point we note that the film susceptibility has been also discussed by Nemirovsky and 
Freed (see Eqs. (3.14d) and (3.16d) in Ref. |§8[). Instead of the (p, 2:1, Z2)-i'epresentation 
of the propagator employed here they used a discrete spectral (p-Kj)-representation. In the 
discrete representation the propagator for Dirichlet boundary conditions is given by 

Gd,M^) = 2 , \ 2 ^ = vr(j + 1)/L, J = 0,1,2,... . (3.34) 
The (p, ^1,2:2)- and (p-/tj)-representation are related by the formula 



2 ^ 

Gd{p, zi, Z2, L,r) = -2_^ sm{KjZi) sm{K jZ2)GD,j{p, r). (3.35) 

j=0 

The one-loop contribution to the total susceptibility is given by 

L L L 

gn + 2 f f f df^-^q 



2 3 y "^"^ J (^ j ^^-^^^ 



Gd{p = 0,zi,z,L,r) GD{q,z,z,L,T) Gd{p = 0, z, Z2, L,t) 
= -f-^— / / f !t Y-i j (t) Y1 sin(K„,2;i) sm{K^^z) 

mi,m2,m3=0 

GD,mAP = t) sin^(K„2z)G'z),m2(g, r) sin(K„3z) fW.{KmzZ2)GD,mi{P = 0, r). 

After performing the integrations one has to evaluate the triple sum. In their calculation 
of the susceptibility Nemirovsky and Freed omitted the terms mi 7^ in the above sum 

20 



which leads to an erroneous expression for the scahng function f{y). If, however, all terms 
in the triple sum are properly taken into account, one obtains, as expected, the same correct 
result for f{y) as via the (p, 2:1, 2;2)-representation. 

The above discussion is focused on the limit y = L/^ 00, i.e., on increasing the film 
thickness at a fixed temperature. In the opposite limit y ^ the film thickness is kept 
fixed and one approaches the bulk critical temperature Tc{L = 00) where ^ diverges as 
C,q{{T — Tc{L = oo))/Tc(L = oo))^'^. For Dirichlet boundary conditions as considered here 
the critical temperature of the film occurs at a lower temperature Tc{L) < Tc{L = 00). 
Therefore the film is not critical at Tc{L = 00) and thus the susceptibility is an analytic 
function of t around t = {T — Tc{L = oo))/Tc(L = 00) = 0. Therefore the finite size scaling 
function g{y — * 0) has the following form: 

g{y ^ 0) = -Ay-''/" - + V + Ey'l'' + (3.37) 



with 



^ = 1^ " + ^^^^ + ^)) + ^^''^ ^^-^^^ 



and 



so that 



/(y ^ 0) = P + Ey^l^ + 0(y2/-). (3.40) 



The numbers 02, &2, and 64 are given in Eq. ( P26|) in Appendix p 3| . For (n, d) = (1, 3) the 



values of the amplitudes to first order in e are V ^ 0.08142 and S ~ —0.01375; for A and B 
see Eq. ( |3.27| ) . The explicit form of the scaling function and its limiting behaviors are given 
in Appendix [B 3| . Figure ^ (a) shows f{y) within mean- field theory and within perturbation 



theory in first order e as well as its corresponding asymptotic behaviors for large and small 
values of y, and Fig. || (b) displays g{y) for large values of y. 
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Our investigations are restricted to temperatures T > Tc. Recently Leite, Sardelich, 
and Coutinho-Filho (LSC) have analyzed amplitude ratios of the specific heat and the 
susceptibility above (T > Tc) and below (T < Tc) the bulk critical temperature in the parallel 
plate geometry for various boundary conditions. These amplitude ratios are functions of the 
scaling variable L/C,± (where C,± is the correlation length above (+) and below (— ) the 
bulk critical temperature) and describe the surface excess and finite-size contributions of 
the system. Their result for the amplitude function of the susceptibility above Tc (see the 
expression for (7+ in Eq. (22) in Ref. [^]) can be expressed in terms of the scaling function 
f{y) as introduced in Eq. ( p. 241) . Within this framework the results of LSC to first order 
in e for Dirichlet boundary conditions are equivalent to the following version of the scaling 
function f{y): 

-1 



fisciy) = y 



-2 



e 

1 - - 
3 



dsfi/2 



y 

s— 

TT 



6yi 



(3.41) 



with 



/i/2(t 



(^2 _ a^yi/2^u 



(3.42) 



exp(27rM) — 1 

For small values of the scaling variable y this scaling function fisciy) deviates even within 
mean-field theory qualitatively from the actual correct form f{y) given in Eqs. ( |B23| ) and 
( p.40| ). Moreover, already for y ^ 10 the difference between Jlsc ^-nd / becomes larger than 
10% in 0{e) and larger then 25% within mean-field theory. These discrepancies are due to 
the fact that even within mean-field theory the results of LSC do not reproduce the correct 
surface excess contributions [p| and finite-size contribution (Eq. (p23|)). 



IV. SCATTERING CROSS SECTION 

A. Scattering theory 

As pointed out in the Introduction the diffuse scattering of X-rays and neutrons under 
grazing incidence allows one to probe the local structure factor near interfaces and in thin 
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films. In this section we discuss how the singularities of the two-point correlation func- 
tion near criticality in a film, as calculated above, translate into singularities of the diffuse 
scattering intensity under the aforementioned experimental conditions. 

We consider a film (0 < 2; < L) composed of a material 2 sandwiched in between 
two halfspaces filled with material 1 (2; < 0) and 3 {z > L), respectively (see Fig. An 
incoming plane wave of X-rays or neutrons with momentum K* = (kj, g^) impinges on the 1-2 
interface at an angle of incidence so that qi = sin a and kj = cosai(cosv5i, sim^j, 0). 
A = 271 / K'' is the wavelength of the X-rays or neutrons. We assume that the media 1 
and 3 are homogeneous and that the 1-2 and 2-3 interfaces are laterally fiat so that their 
contributions to diffuse scattering can be ignored. Within the plane of incidence there is 
a specularly reflected wave with K'' = (kj,— g^). The mean value of the electron density 
in the case of X-rays and of the scattering length in the case of neutrons determine the 
intensity of the reflected beam whereas fluctuations around the mean value give rise to 
scattered intensity in off-specular directions K-'^ = (kj,gj < 0) with gj = — i^'-'^ sin and 
kj = K-^ cosa/(cos v?/, sin 0). We consider only elastic scattering, i.e., = = = 
K. (For the more complex case of neutron scattering under grazing incidence from magnetic 
systems see Ref. P5|.) 

In order to proceed we assume that the mean values of the electron density or of the 
scattering length density in each medium is constant and varies steplike across the two 
interfaces 1-2 and 2-3. This gives rise to the following indices of refraction |50| : 



^<0:n = ni = l (4.1) 
Q<z<L: n = n2 = l — 82 + «/?2 
z > L : n = = I - 5^ + i/S^. 

In Eq. ( [4.1| ) we consider the case that medium 1 is vacuum and the generic case for hard 
X-rays that Re < 1 in condensed matter. Although for neutrons one can also have 
Re n > 1, in order to limit the number of possible relative values of the indices of refraction 
for the materials 1, 2, and 3 we do not analyze this latter case in more detail. For X- 
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rays 5 = ^"^^Yli^i^i the extinction coefficient /3 = ■^Yli^i^a,i = wfiere re = 
. 2 = 2.814 • 10"^ A is tfie classical electron radius, iVj the number density of atoms of 
species i with Zj electrons and absorption cross section aa,i- For neutrons 5 = ^ Ylii^i^i 
and P — X^j -/VjiTt^i where 6j is the nuclear scattering length of species i. at,i is the cross 
section taking into account incoherent scattering and nuclear reactions. Typically 5 and /3 
are of the order 10~^. For Re n <1 total external reflection occurs for a < etc- For L = 00 
one has — (252) whereas for L = ctcia — (2^3)^/^. Since the angle of total reflection 
depends only on the difference n{z —00) —n{z +00) > 0, for any finite < L < 00 the 
incoming wave is totally reflected for a < ctcia, independent of the index of refraction within 
the film. Nonetheless the types of waves propagating in the film depend on whether a ^ aci2 
(see below). For the present setup the wave field has the form ^(r, K*) = e*^^ '"11^(2;, a) with 

jjjiz, a) = < s+(Q;)e^«2(")^ + s-{a)e-'i^^''> , 0<z<L (4.2) 
tL(a)e'«3(")^ , z> L 

where 

TLia) = ((gi - q2){q2 + Qs) + e^'^'^Qi + q2){q2 - q3))/A{a), (4.3) 
s+(a) = 2qi(q2 + q3)/Ha), 
s_{a) = 2q,{q2-qs)e''^^'^/A{a), 

ti(a)=4gi52e'(«^-^^)VA(«), 

A(«) = {qi + q2){q2 + gs) + e2^^^^(gi - 52) (^2 - qs). 

Since the scattering cross section is independent of the intensity of the incoming beam 
without loss of generality we have set the amplitude of ^'(r,K*) equal to 1. The vertical 
components of the momentum are given by 

qi{a) — Ksina, (4.4) 
qj{a) — K^n'j — cos^ a ~ KyJ sin^ a — 25 j + 2i(5j 



— K\l sin a — sin ctcij + 2ij3j , j = 2, 3. 
24 



In the limiting case that the film turns into a semi- infinite substrate, i.e., L = oo one has 

(4.5) 

too/2 (a)e*^2(")- , 2 > 

with 

^oo/2(a) = {qi - q2)l {qi + 92), (4.6) 
too/2 (a) = 2gi/(gi + 92)- 

The vertical momentum components qj{<y) have a positive imaginary part which is due to 
the extinction coefficient l3j for a > adj and which is present for a < Odj even in the 
absence of absorption. This gives rise to an exponentially damped evanescent wave with 
a penetration depth Ij = {Im qj{a))^^ which increases steeply for a y adj and would 
diverge if Pj = 0. Within the film there is a superposition of two fields s+{a)e^'^^^°'^^ and 
s_(a)e~*^2*^°^^ (Eq. (|4.2| )); in the three cases a < (Xd2 and = 0, a > ad2 and 7^ 0, 
and a < ad2 and (32 7^ 0, q2{(y) has a nonzero imaginary part leading to an exponentially 
increasing and decreasing contribution for increasing z. The decreasing part corresponds to 
the damping of the incident wave whereas the increasing part corresponds to the damping 
of the refiected wave generated by the interface 2-3. 

Equation (|4.2| ) describes the wave field \I'(r, K*) in the absence of any fiuctuations. This 
wave field is scattered at the fiuctuating inhomogeneities within the film giving rise to diffuse 
scattering intensities in off-specular directions. 

The computation of this intensity requires one to specify the nature of ffuctuations. In 
the present context this amounts to specifying the kind of system undergoing the continuous 
phase transition in the film and to choose the appropriate order parameter. As described 
in the Introduction the most promising candidates for these kind of phenomena are binary 
alloys undergoing a continuous order- disorder phase transition concerning the occupation 
of fixed lattice sites {R^}. (Magnetic films are equally well suited. However, the magnetic 



scattering of neutrons P3 or of X-rays is more complicated and requires separate analy- 



ses. Although the details will differ from the analysis given below, the key features of the 
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singularities are expected to be born out similarly.) In these systems a given configuration 
is characterized by spin-type variables {5*; = ±1} such that 5*^ = +1(— 1) states that the 
lattice site R/ is occupied by a B{A) atom. Accordingly the number density of electrons for 
such a configuration is 

P(r) = IJ2 {^^(^ - + ^^(^ - + [/^^(^ - - ^^(^ - (4.7) 

I 

where pa{b){y) is the electron number density in a single unit cell Vceii occupied by an A{B) 
atom. (In the case of neutron scattering p(r) stands for the scattering length density and 
Pa(b){t') = &A(B)'^(i") where is the mean scattering length of the nuclei of species A{B).) 
The ordered state of this system corresponds to a configuration in which the sign of Si alter- 
nates from one lattice site to any of the neighboring ones. In this ground state the staggered 
"magnetization" OPi = Sie^'^"''^'- is spatially constant if the reciprocal lattice vector of 
the sublattice structure is chosen such that e*Tm-(R-i-R-;') = _i fQ^- nearest neighbor sites 
R;,R/'. In the reciprocal space the positions of the reciprocal sublattice vectors are 
halfway in between the reciprocal lattice vectors Gm with e**^™'-^' characterizing the un- 
derlying lattice structure of the solid. (For the sake of simplicity as far as the scattering 
theory is concerned we do not consider here explicitly the case of systems like Fe^Al whose 
description requires the introduction of several sub lattices.) Upon approaching the critical 
temperature of the continuous order-disorder transition the thermal average {OPi) vanishes 
qualifying OPi as an appropriate order parameter. 

In the critical contribution to the bulk scattering cross section a nonzero value of (OPi) 
leads to superlattice Bragg peaks |^ : 

V / cell ^ 

where is the classical electron radius, K.^/K are the directions of observation, e is 
the polarization vector of the incoming electromagnetic wave, F = {Fa — Fb)/2 where 
Fa(b)(K) = Jy d^r p A{B){^)e'^^'^ is the atomic form factor of the atom A{B), e"'^^^^ is the 
Debye Waller factor and Ny is the number of lattice sites in the sample. With the inde- 
pendent knowledge of all prefactors in Eq. ( [4.8|) the asymptotic temperature dependence 
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Bragg 



e 



bulk 



of yields (OPi) = B'{—tY. As discussed in Sec. |T| this experimental value for 



B' enters into Eq. ( p.21|) and there replaces B if G{p, Zi, Z2, L,t) corresponds to the pair 
correlation function {OPiOPi') as considered below. Similarly the singular diffuse scattering 
around a superlattice Bragg peak is given by 

diffuse /-TT-t \ 2 



dn) 



bulk 



/TZ-f \ _ |2 

( X ej Fe-^l J2 ((OPiOPi') - {OPi){OP,)) e^^^^^-""^'^ (4.9) 

Gbulk{<i, t) 



K 



X e 



Fe 



R;,Rj/ 

2 Ny 



cell 



with q = — K' — Tm- In the second part of Eq. ( |4.9| ) we have performed the continuum 
limit replacing the lattice sums by integrals (see Eq. ( |A3| ) and the last paragraph in Sec. ^ 
because for ^ ^ 00 the lattice structure becomes irrelevant. From studying the temperature 
dependence of Eq. ( [4.9|) for T > one can infer the correlation length C, and its amplitude 
introduced in Sec. ||. We note that for q small compared with the inverse lattice spacing 
a Eqs. ( |A^ ) and (|A5|) can be applied to Eq. ( [4. 91) provided B is replaced by B' as determined 
from Eq. ( ^ . 

Equipped with this knowledge about the critical bulk scattering (i.e., above the angle 
of total reflection and for a bulk sample) we can now turn to the critical diffuse scattering 
from the film. Within the so-called distorted wave Born approximation and for the model 
of the film as described above one finds the following expression for the singular part of the 
coherent scattering cross section ||49|| : 



da 

dn 



(4.10) 



^W^ceii illuminated surface area where A^y is the number of lattice sites at 



dzi J dZ2lljf{Zi)llJi{zi)llj*{z2)lp}{z2)G{p,Zi,Z2,L,t) 


where A\ 

the surface and V^^n is the two-dimensional unit cell of the surface, a is the lattice spacing of 
the cubic lattice, ipij^z) = ipi^z, a = aij) (see Eq. ( [4.2| ) and Fig. |T0|), and p = k-^ — k* — Tm 
assuming that the film surfaces are cut such that Tm is parallel to them. G is the lateral 
Fourier transform of the two-point order parameter correlation function: 
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G{j>,z,,zM)=^-f^ e^P-(^ii -^11 ) (4.11) 



^(m) ^{m') 



X ((OP(r[|'"\zOOP(r[|'"'\z2)) - (OP(r(|"^\ zi))(OP(r;|'"'\ ^2)) 
y e-'P"ii G(r||,2i,Z2,I^,t) 



on the lattice and in the continuum hmit, respectively. Thus after replacing the nonuniversal 
amplitude B in Eq. ( |2.21|) by B' as obtained from Eq. (|4.8|) for (QP^) we can study the 



scattering cross section in Eq. (|4.10| ) by using all the information about G{p, Zi, Z2, L,t) 
obtained in the previous section, provided all lengths and 1/p are sufficiently large compared 
with the lattice spacing a so that the continuum description is applicable. 

In view of the properties of the wave functions ip {~ only their functional forms for 
< z < L enter into S (see Eq. (^I^))— ) and of the scaling form for G{p, zi, Z2, L,t) (see 



Eq. ( |2.21| )) one has for a^j, aci2,ci3 ^ 1 and j32,3 = 0: 



L li If ai aci2. 



E = B'i^^r^'n { ^ ) a{p^, J, J, ^f, ^) (4.12) 



where the dimensionless function a is given by (Eq. ( |2.21| )) 
1 1 



a = y dxi J dx2ipf{zi = xiL)ijji{zi = xiL)ip*{z2 = X2L)^jjj{z2 = X2L) (4.13) 


xl~'^gnipLx2, jXi, ^X2, X2). 

The two variables p^ and L/$, of cr stem from the scaling function of the pair correlation 
function whereas the dependences of a on k/L, If/L, ai/aci2, and acu/c^ds are due to the 
wave functions. For aij < acvi 

/(2) 

k< = I ° (4.14) 

correspond to the penetration depths of the incoming (z) and outgoing (/) evanescent wave, 
respectively, within the film material 2. 1^^^ = {Kacu)'^ is the minimal penetration depth 
lij{aij = 0) in the film material. Typically Iq is of the order of 30 A [Q. For > aci2 
and af > aci2 the corresponding quantities k and respectively, are purely imaginary. 
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B. Interplay of length scales 



The scattering cross section reflects the rich interplay of flve length scales: ^, li, 
If, and L. Scahng reduces that to four independent scaling variables; moreover there is a 
parametric dependence on ai/acu and on the material constant aci2/(^ci3- It is beyond the 
scope of the present analysis to provide an exhaustive discussion of the full dependence on 
all these variables. Instead we discuss some general aspects and analyze a few specific cases 
in more detail in order to highlight the key features of the diffuse scattering intensity. The 
following cases have to be distinguished (for T >Tc): 

I a) kj <^ L and total reflection at 1-3 interface: ^ is proportional to the scattering 
volume A\^'min{li, If) 

1. ^ <^ hj <^ L: bulk behavior convoluted with evanescent waves 

2. ^ ~ lij <^ L: crossover bulk / y surface behavior convoluted with evanescent 
waves 

3. lij <S ^ ^ L: Y surface behavior convoluted with evanescent waves 

4. lij <C ^ ~ I/: Y surface behavior plus distant wall correction convoluted with 
evanescent waves 

5. lij <^ L <^ ^: fllm behavior near one wall convoluted with evanescent waves 

I b) kj <^ L and no total reflection at 1-3 interface: the difference to I a) is exponentially 
small, i.e., ~ e~^^ . (The volume contribution to ^ from material 3 is insigniflcant 
because it does not exhibit critical fluctuations.) 

II a) lij ~ L and total reflection at 1-3 interface: crossover between ^ ~ A\\min{li,lf) to 




1- e < kf 



L: Y surface behavior convoluted with fllm wave functions 



2- e ~ kf 



L: crossover bulk / ^ surface behavior convoluted with film wave 



functions 
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3. lij ~ L ^ ^ — s> oo: film behavior convoluted with film wave functions 

II b) lij ~ L and no total reflection at 1-3 interface: crossover ^ ~ A\\min{li,lf) A\\L 
(Again, the volume contribution from material 3 is regarded to be insignificant and is 
not taken into account.) 

Ill a) lij ^ L and total reflection at 1-3 interface: ^ ~ AyL 

1. ^ L ^ lij: bulk behavior convoluted with film wave functions 

2. ^ ~ L ^ lij: crossover between bulk and film behavior (including two surface 
contributions and distant wall corrections) convoluted with film wave functions 

3. L ^ ^ ^ lij: film behavior convoluted with film wave functions 

4. L <^ ^ lij: film behavior convoluted with film wave functions 

5. L <^ lij <^ ^ —>■ oo: film behavior convoluted with film wave functions 

III b) lij ^ L and no total reflection at 1-3 interface: ^ ~ AyL (The volume contribution 

from material 3 is regarded to be insignificant.) 

IV a) lij imaginary and total reflection at 1-3 interface: ^ ~ AyL 

1. ^ <^ L: three-dimensional bulk behavior probed by undistorted plane waves 

2. L ^: film behavior probed by undistorted plane waves 

IV b) lij imaginary and no total reflection at 1-3 interface: ^ ~ AyL (in addition to an 
insignificant volume contribution from material 3). 

C. Susceptibility from the scattering cross section 

For large penetration depths lij ^ L the product of wave fields in Eq. ( [4.13| ) is approxi- 
mately constant. In this case for p = the universal scaling function a of Eq. ( [4.13| ) reduces 
up to a prefactor to the scaling function / of the total susceptibility (see Eqs. ( ^.24| ) and 
(PI)), i.e., 
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a(^)=.(p^ = 0.^.^i = oo.^ = oo.^<l.^) ~/(^). (4,15) 

\U \ i L L aci2 aci3/ \i J 

In this limit the dependences on ai/aci2 and on aci2/«ci3 drop out for < Ods; for a, > 
there is an insignificant bulk contribution from material 3. The five different cases 1. — 5. in 
III a) are characterized by the various contributions of asymptotic behaviors to the scaling 
function a{y = L/0 ~ fiv) (see Eqs. and (lOOl) ), i.e., bulk: Ay^^^"^, surface: By-^+'', 

distant wall: Cy~^^'^e~y , and film behavior: T) + Sy^^'^ . In Fig. |TT] we show the normalized 
scaling function of the scattering cross section a^iy) = cr(?/)/cr(0) (Eq. ( |4.13| )) within mean- 
field and within first-order perturbation theory as well as the asymptotic behaviors of the 
normalized scaling function fo{y) = /(y)//(0) of the total susceptibility f{y — * 0) (Eq. 
( p.40| )) and f{y ^ oo) (Eq. ( p.26| )) using mean-field exponents and amplitudes and best 
values for the exponents and amplitudes to first order in e, respectively. The cases III a) or 
b) with lateral momentum p = are the appropriate scattering setups in order to measure 
the various asymptotic behaviors of the total susceptibility by varying the temperature. 

Figure |l^ shows the ratio of the normalized scaling functions fo{y)/cro{y) = ^^j^ 
all four cases I a) — IV a) within mean-field theory (MET) and within perturbation theory 
(PT), respectively. For large penetration depths kj ^ L (see case III a)) the deviation of 
the scaling function ao{y) of the scattering cross section from the scaling function fo{y) of 
the total susceptibility is small (see the solid lines in Fig. ^ (a) and (b)). If the penetration 
depths are of the order of the film thickness, lij ~ L (see case II a)), the wavefields in Eq. 
( [4.13| ) contribute and the deviation from the total susceptibility becomes visible at large 
values of the scaling variable y. For y ^ and y ^ oo the dotted lines attain constant 
values so that there are the same critical exponents but different amplitudes for the leading 
asymptotic behaviors of do and /q. If the penetration depths are smaller than the film 
thickness, kj <^ L (see case la)), this deviation is much more pronounced (see dashed 
lines). The difference in the amplitudes is decreased if aij > acu, i-e., for imaginary Zjj 
(see case IV a) and dashed-dotted lines). 
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D. Dependence on the film thickness 



In order to reveal the (1/i^) ^'^''n cusp singularity in the scattering cross section. We 
consider the case p = t = and introduce the corresponding scattering function 

L L 

^L = j dzi J dz2^f{zi)^i{zi)ij*{z2)rf{z2)Gip = 0,Zi,Z2,L,t = 0), (4.16) 





where the wave fields are given in Eq. ( [4.2|) . For the correlation function G we use the 
asymptotic expansion given by Eq. (|3.21| ). Furthermore we introduce the scattering function 
of the semi-infinite system 



oo oo 



= I dzi dz2%l)^iJ izA^j^zMT {z2Wl/^ (4.17) 

2 / / 2 2 2 2 





with the wave fields and the correlation function given in Eq. (|3) and Eq. ( 3.18|) , respec 
tively. The ratio of Eqs. ( 4.16 ) and ( 4.17| ) defines the scattering function 

Si 



(4.18) 



for p = t = 0, where the film thickness L and the momentum K of the scattered wave form 
the scaling variable, the angles a = aci2,ci3} characterize the scattering geometry, and 
the extinction coefficients (3 = {(32, (3^} take into account photo absorption. From Eqs. (|C3|) 
- (|C8|) in Appendix one obtains the asymptotic expansion 



/ 1 \ 

S{LK 00; a, P) = so{LK; a, p) + Si{LK- «, /5)Ci (^— j + . . . 



(4.19) 



with 



So{LK ^ OO] a, (3) ~ 1 + 4^^(Lir; a, /3)e-^^"o"("'^), 
si{LK 00; a, (3) ~ sf\LK = 00; a, (3) + s^P{LK- a, f3)e~^^'''?^'''^\ 



(4.20) 



and Ci given by Eq. (|3.9|). The functions sq and si carry the L-dependence of the wave 
functions (see Appendix 0). The L-dependence due to the correlation function is given by 
the cusp singularity Ci{l/ LK)~^^^\\ . The range of the values of the scaling variable {LK)~^ 
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is limited by the validity of the continuum theory applied here, i.e., L > 30 A and the 
distorted-wave Born approximation, i.e., K > 1 A^^, leading to {LK)~^ < ^. For small 
angles, i.e., for grazing incidence scattering experiments Eq. ([4.4|) reduces 



qi{a) ~ Ka, (4.21) 



qj{a) ~ - + 2i/3j, j = 2, 3. 

Photo absorption, 7^ 0, or evanescent scattering, ctjj < acu, turn q2 into an imaginary 
quantity, which leads to a real part of Sq^^ and s^i^ in Eq. ( [4.20| ). If at least one angle ai or a f 
is larger than the critical angle the functions s^q^ and sf'^ have a real and an imaginary 
part. In the latter case one expects that the scattering function S in Eq. ([4.19|) exhibits 



an oscillatory behavior. In Fig. we show the exponentiated scattering function and its 
asymptotic form for various scattering geometries. The exponentiated form is obtained by 
subtracting the leading behavior of the one-loop e-expanded result (defined by Eq. ( |4.18| )) 



and by adding the leading behavior (see Eq. (|C^)) calculated with the best available critical 
exponents (77 ~ 0.031, t]± ~ 0.75, and r^y ~ 1.48). The dashed line in Fig. |13| corresponds 
to the leading asymptotic behavior, if the L dependence of the wave fields is neglected: 

S{LK ^ 00; a, P) = 1 + s^°\lK = 00; a, P)Ci(-^) + . . . . (4.22) 



Thus the full lines in Fig. |I3| take into account the whole L dependence stemming from both 
the scattering theory and the correlation function, whereas the dashed lines take into ac- 
count only the leading asymptotic L dependence of the correlation function. The oscillatory 
behavior appearing for at ^ acu ^ otf stems from the scattering theory (see Fig. p!3|). 



For the case aij < aci2,ci3 in Fig. |13| half of the maximum value of the scattering 
function S is reached for LK ~ 1.5 ■ 10~'^. This corresponds to a film thickness L ~ 600 A, 
i.e., 200 monolayers (with K c^l A^^ and 1 monolayer ^ 3 A thick); 90% of the maximum 
value of 5* is reached for LK ^ 5 ■ 10^^ which corresponds to a film thickness L ~ 20000 
A or 6700 monolayers. This demonstrates the slow convergence to the semi-infinite limit. 
The spatial resolution is determined by the uncertainty of the film thickness. With AL ~ 3 
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A (1 monolayer) this gives KAL ~ 3 leading to a resolution of A{LK)~^ ~ 3/{LKy, which 
is not visible on the scale of Fig. |13[ Based on these considerations we conclude that the 
oscillations are experimentally accessible. 



E. Emergence of cusp singularities 

In the following we analyze how the cusp singularities emerge in the limit of vanishing 
scaling variables. To his end we chose as an example a scattering function of the two scaling 
variables p/K and LK. Analogous to Eq. (|4.16| ) we define the quantity 

L L 

T.p^L = j dzi j dz2i'f{zi)iJi{zi)iJj*{z2)tp}{z2)G{p,zi,Z2,L,t = 0). (4.23) 



Together with Eq. ( [4. 17] ) this leads to the scattering function 



S{p/K,LK;a,(3) = ^ (4.24) 



2 



where a = {aij, aci2,ci3} denotes the set of angles and /3 = {^2,3} the extinction coefficients. 
As in Subsec. [Ill A| and Eq. ( p.l3| ) we introduce polar coordinates 



Lu = \/{pIKY + (Li^')-2, ^ = arctan(pL) , (4.25) 
(LK)~^ = uj cos if, p/K = uj simp . 

This leads to the relation 

S{p/K,LK;a,(3) = 5'(a;sinv3, {u cos (p)~^ ; a , (3) = Spoiar{u},(p]a, (3) (4.26) 

so that the leading asymptotic behavior is given by 

SpoiariuJ ^ 0, ip; «, /3) = Soi^p; «, /3) + S.icp; a, /3) + . . . (4.27) 

with So{ip; a,P) = 1. The amplitude 5*1 of the leading asymptotic behavior cj^^+'^ii depends 
not only on the polar variable (p, as it is the case for the corresponding correlation function 
(see Subsec. [Ill A| ), but also on the parameters a and /? characterizing the scattering pro- 



cess. Within mean-field theory this amplitude is defined in Appendix [C 2| by Eq. P19[ . In 
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Fig. ^ (a) we show the exponentiated scattering function S{p/ K, LK; a, P) (Eq. ([4.24|) ), 
where we have subtracted the leading asymptotic behavior from the mean-field expression 
of the scattering function and added the exponentiated form. (See Eqs. (|C11| ) and ( pi9|) 
in Subappendix |C 2| ) ; the scattering function S (Eq. ( |4.24| )) is a sum (see Eq. ( pij ) in Ap- 



pendix ^ of functions of the type S as discussed in Eq. ( |(J11| ) in Subappendix |C 2| .) Figure 
[14| (b) illustrates the emergence of the {p/ K)~^^'^^^ cusp for increasing film thickness, i.e., 
(LK)^^ 0. Figure Fig. [l^ (c) shows the emergence of the {1/{LK))~'^~^'^\\ cusp for vanish- 
ing lateral momentum p/K 0. In the later case the vertical cross sections of the manifold 
are not monotonous; they exhibit a maximum (•) at 1/L ^ 0. Figure |l^ corresponds to 
scattering angles a^j < aci2,ci3 which yields a monotonous behavior of the scattering func- 
tion. Analogous considerations describe the emergence of the cusp singularities in the ^-L 
and ^-p dependences (see Subappendix 



V. SUMMARY 

By using fieldtheoretic renormalization group theory we have studied the singular part 
of the two-point correlation function in a film of thickness L near the critical point of the 
corresponding bulk system. For T > and Dirichlet boundary conditions we have obtained 
the following main results: 

(1) The two-point correlation function as a function of the lateral momentum p corre- 
sponding to the d — 1 translationally invariant directions of the film geometry, the coordi- 
nates Zi and Z2 perpendicular to the parallel surfaces, the film thickness L, and temperature 
t = (T — Tc)/Tc (or equivalently the bulk correlation length ^ = C.Qt^'^) exhibits three cusp 
singularities: ^^-'^+'''11 for t = and L = oo, (l/^)~^+''ii for p = and L = oo, and (l/L)""'^^'''!! 
ioT p = t = (see Eqs. ( |3.19| ) - ( |3.21| ) and Fig. |T]). The emergence of these three cusp 
singularities is revealed by studying appropriate scaling functions of two scaling variables 
(see Eqs. ( |330| ) - (|32|) and Figs. |- |). 



(2) The film correlation function calculated up to first order perturbation theory in 
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e = A — d satisfies the so-called product rule derived by Parry and Swain for the correlation 
function algebra of inhomogeneous fluids in Ref. (see Eq. ( |3.22| )). 

(3) By setting p = and integrating over the perpendicular coordinates zi and Z2 we 
obtain the total susceptibility of the film (Eq. ( |3.23| )). Its dependence on L and ^ is described 
by a universal scaling function f{y = L/^) (see Eq. ( |3.24| ) and Fig. |^) and exhibits a typical 
film behavior: f{y) is analytic for y ^ and f{y oo) contains the bulk-, surface-, 
and finite-size contributions (see Eqs. (|3.4CI|) and (|3.26| ), respectively). These properties 
are similar to those of the specific heat of a critical film Our results correct previous 
findings in the literature |]68| , |82| (see the discussions of Eqs. (|3.36|) and ( p.41|) ). 

(4) In view of proposed experimental tests with X-rays and neutrons under grazing 
incidence (see Fig. |10|), as discussed in detail in Sec. |, we have calculated the critical 
diffuse scattering from the film within the so-called distorted wave Born approximation . 
The scattering intensity is a function of the lateral momentum transfer p, film thickness L, 
bulk correlation length ^, penetration depths lij of the incoming (z) and outgoing (/) waves, 
the critical angles of total refiection and Ocis and the extinction coefficients and /^a 
of the film (2) and of the underlying substrate (3) (see Fig. pUj). 

(5) For various ratios of L, ^, and lij the scattering function shows the crossover between 
analytic, bulk, surface, and finite-size behavior (see Figs. |ll] and |12|). By varying the tem- 
perature, a scattering experiment for p = and lij ^ L gives access to the aforementioned 
scaling function f{y) of the total susceptibility (Eq. ( [4.15|) ). 

(6) For p = t = the leading singular behavior of the scattering function is given by 
the cusp singularity (l/LK)^^^''ii, where K is the momentum of the incoming wave (Eq. 
( [4.19| )). The maximal scattering intensity for L — > oo is reached only very slowly. For certain 



scattering geometries the L-dependence exhibits an oscillatory behavior (see Fig. |T3]). 

(7) The film thickness and momentum cusp singularities of the correlation function are 



borne out in the scattering cross section and are analyzed in Fig. 14 
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APPENDIX A: AMPLITUDES 



The amplitudes of the singular behavior of bulk correlation functions are nonuniversal. 
There are two independent ones in the sense that any two of them allow one to express 
any other in terms of these two and universal amplitude ratios [79,^|. As one of these 



nonuniversal amplitudes in Sec. || we have introduced and fixed the amplitude of the 
bulk correlation length (see Eq. (|2.8| )). Other nonuniversal amplitudes are given by the 
temperature dependence of the mean value of the field 0(x) below Tc, 

(0(x)) = 5(-t)^ (Al) 

by the decay of the two-point correlation function in real space at for large distances 
|x — x'l, 

(0(x)0(x')) = /^|x - x'|-('^-2+^), (A2) 
and in momentum space for small g, 

j d'^a;e*'i (''-^')(0(x)0(x')) = Gu.ik{q,t = 0) = bq-^+\ (A3) 



where 



D/D = X = 22-%^/2^lL^Z/^L^. (A4) 
/ r(rf/2 - (1 - r//2)) ^ ^ 



D can be expressed in terms of B, ^q, and a universal number R [ [r^ , pO| 



D = B^^^y-^+m (A5) 
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with R = RcQs/iRp^. For {n,d) = (1,3) one has Rc ~ 0.066, Q3 ^ 0.922, and R+ ~ 0.27 
|7|,|80| leading to R - 3.09. 

A Fourier transformation in the z-direction of the bulk correlation function Gbuik{(l) = 
jjq-'i+n with = + k'^ is leading to its p-2;-representation 

/°° dk ^ g«fe(2i--z2) 

f) /-oo iKp{zi-Z2) 

For — 2:2) ^ this leads to 



^^"''^^^^-^ 27^ r(l-ry/2) • 



(A7) 



In the limits L — » 00, zi + Z2 ^ 00, ^ — 00, and p{zi — Z2) the two-point correlation 
function in the film reduces to its bulk form. According to Eqs. ( |2.12| ) and ( p.l7| ) this 
implies 

D r(l/2-V2) _ „2.,+^.-2+, R r(i/2-V2) 
- 20F r(l - V2) - ^ 20F r(i - V2) ^""^^ 

For the three-dimensional Ising model the universal number U has the value U ~ 1.58. 

The knowledge of the perturbative result for G{p, zi, Z2, L, t) (see Subappendix |B 2|) en- 
ables one to express the nonuniversal amplitudes Qx, x = 1 - IV, in terms of Qy- For exam- 
ple the universal ratio Qu/Qv is determined by the normalizations of the scaling functions, 
i.e., c/„(0, 0,0,0) = 1 (Eq. (gJ^)) and ^v(oo, 0, 00, 00) = 1 (Eq. (|2T7|) ). The e-expansion of 
this ratio is given by 

^„/6?v = 2(1 + + Oie')). (A9) 

V n + 8 / 

The amplitudes Q^, x = I - IV, can have bulk, half space or film character, depending 
on the normalization limits of the scaling functions Qx- Gn and Q^, are half space and bulk 
amplitudes, respectively. Bulk amplitudes are independent of the boundary conditions, 
halfspace amplitudes depend on the boundary condition of the surface, and film amplitudes 
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depend on the boundary conditions of both surfaces. Combining Eqs. 
we arrive at 







1 + e 



n + 8 ' V"0F r(l -77/2) 
With R ~ 3.09 (Eq. (^)) and r] ~ 0.031 one has for the 3d Ising model 



6;„ = 7^E2(eo+)'^-2+^ 



8D, and 



(AlO) 



(All) 



APPENDIX B: ONE-LOOP RESULTS 



1. Correlation functions for zi = Z2 



With the abbreviation e = e^^, so that e = |, |, ^ for the Ising, XY, Heisenberg model 
in = 3, the renormalized two-point correlation function in one- loop order (Eqs. (|2.3|) - (|2.6| )) 



is given explicitly as (see also Ref. ||84 



G{p, z,L = 00, t = 0) = Guz'-^giiu = pz) (Bl) 

/I _ -2u ~ \ 

= ( + ^ [-1Ei{-1u) + e'"Ez(-2M) + e-'"E^(2M)) + 0{^) \ . 

Ei{x) is the exponential integral function. In accordance with the normalization gi{0) = 1 
this yields Qn = /i~''(l + e + 0{e^)). 

The temperature dependence is described by the scaling function g2{v) with (72(0) = 1: 



Gip = 0, L = 00, t) = gnz'-''g2iv = z/0 (B2) 
-.-.-^(^p + |((e--l)iro(2.) 

^ to {m2k+i) j + ^^^V- 

\l'(x) and -^'0(3;) denote the psi and Bessel function, respectively |^ (see also Ref. [^ ). 
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Finally, the dependence of the critical structure factor on the film thickness is governed 
by a third scaling function g3{w), ^'3(0) = 1, < w < 1: 



G{p = 0,z,L,t = 0)= Gnz^-'^gsiw = z/L) 



(B3) 



^1 - w + e| - -^'^O- ~ '^Y 

— (1 — 2w) l + C£; + lnwH ■ 

V w 

+ (1 - W7) (2 + Ce + Inw; - ^2tiH - A^H) } + ^(e')) 



with the abbreviations 



Sk,Aw) = 2^ ^ . h yw) 

n=k 



dx e ^'^ lb 



(B4) 



n\{n -I) ' 7i - 1 

Bn{w) are Bernoulli polynomials For the critical structure factor in the semi-infinite 

system one has 



G{p = 0,z,L = oo,t = 0) = g,,z^~" = fi-'^z^-" (l + e + 0{e^ 



(B5) 



From the explicit forms for Qi, i = 1,2, 3, in Eqs. (|B1D - ( P3|) together with 7]\\ = 2 — e + (9(e^) 
one infers the following limiting behaviors: 



u 



l + e(l-Cs-ln2) + C(e^ 



(B6) 



and 



g^{v ^ 0) = 1 + + Oi: 



5i 



l + l{l-CE) + 0{t^) 



(B7) 



g^[w ^ 0) = 1 + CiW-^+^ll + 0{W^) 



2 

1 + ~Ce + 2(^2 + h) - 1) + 0{e' 



(B8) 
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where Ce ~ 0.5772 is Euler's constant. 5*2 is given by a sum over Bernoulli numbers and /i 
by an integral: 



Br, 



n=2 



n\{n — 1) 



~ 8.2877 -lO-^ Ji 



1 1 

dx ~ 0.2868. 

— 1 X 



(B9) 



For the exponentiation of the scaling functions hi{u,v), h2{u,w), and h^{v,w) we have 
calculated the amplitude functions h\^\lp)^ H[^\(p), and H[^\ip) (see Eqs. (|3.14|) and 
( p.lSp ). Their e-expansions are 

Hi'\^) = 

ip = arctan((p,^)~^) 



^, , sin 09 cos 09, 1 + cos 09 \ 2\ 

In— 11 + —^ In- ^ + a,) + 0{e') 

2 2 1 — cos (f / 



(BIO) 



1 + e 



2 tan ip 



^ g2 tan ip 



1 + e 



2 tan 1/3 



]^ g2tan(/3 



^1 — Ce — ln(2 sin ip) — To{(p) 



cos if 



+Ti(09) cot 09 H COt^09) ; 

^ 12 V 3(1 _ e2tan</,)(^l _ g-2tan^) 

ip = arctan(pL), 



(Bll) 



and 



_j_ g2tan93 
1 g2tan(/p 

siny9 



_/ 1 _|_ g2tan93 , 

+ g(sin(y) ^ _^^^^^^ (^1 -Ce- ln(sin (^) + J+((^) + Ji-(v9) 

(27r + 8J°(<^) tanip)) + 0{^) , (B12) 



1 g2 tan 99^(1 g— 2 tan 1/3 

99 = arctan(L/^), 



with ai ~ 0.2704 and the integrals 



2^0 (V') 



1 1 

+ 



g e* — 1 Vt + 2tan<y9 ' t — 2tany9 
dt ( t t 



Q e* — 1 Vt — 2tany9 t + 2tan99 



dt^rn • 

g2i tan ip 1 ' 



and 



1 g2nan^_l^±l 



(B13) 
(B14) 
(B15) 

(B16) 
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2. Correlation function for z\ ^ 



This is the most general case from which all results given above can be derived. We 
present -Zi, 2^2, ^) in terms of the scahng function (Eq. 



G{p,zi,Z2,L,t) = GjP ^''gi{x = p^,u = zi/^,v = Z2/^,y = L/Cj 



(B17) 



rp ( p-a{u~v) I p—a{v-u) „-a(ji+j;) „a(ji+'u) 



2a r ' e2y"-l 

+e( Jo(a;,u,f,l/) + Jn{x,u,v,y) + Ji{x,u,v,y)) + C(e^) 



with 



^^^^?73t{^""'""''(1 + "I" - ^1) - + + ^)) (B18) 

( e-"("-^)(l + a(« -v) + -^^) + e-'^(''-")(l + a{v - u) + ^^"^ ^ 



1 _ e-2s/a 

2ya 



2ya 



J^{x,u,v,y) 



TC X 



(B19) 



and 



J,{x,u,v,y) = -l^^e--^^-^^)^:^^ + e'^(^+^-^^) 



1 X 
4a"2 



1 _ e~^^" 



(B20) 



with a = ^/T+x^ and 



-2xis I — 2x2S 



1 1 



s — a s 



rg-2{y-a;i)s _|_ g-2(?/-X2)s' 



(B21) 
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3. The susceptibility 



The one-loop result of the total susceptibility (Eq. (|3.23| )) for Dirichlet boundary condi- 
tions is given by (7s = 7 + t^) 



x{L,t) = B\aY+^n(^Y''' f{y = L/0 



with 



(B22) 



f{y) = y-' 
+ 



4 1 f4 1 2 / 1 + e-y 

ei — — ^ + -— —n{ 1 



yey + 1 iyey + 1 y{l + e-yy 



+ 4 + 8- 



12 1 - e-y 



1 + e-yy yl + e-yJJ e^^y - I 

1 



ds 



+ - 



2 1 - 



ds 



1 



1 2 

+ 



yl + e~y J e'^'y-lXs-l s + 1 s + 1/2 s - 1/2 



(B23) 



In the limit y ^ oo the two integrals entering into Eq. (P23|) vanish ~ e~^^ and therefore 
they do not contribute to the terms considered in Eq. (|3.26| ). However, in the limit y 
these two integrals contribute to the terms considered in Eq. ( p.40|) : 



= - -^y + a2 - - Iny + Oiy) 



(B24) 



and 



with 



My) 



°° , Vs^^/ 112 2 

ds— \ ; — 

1 e^'y-l\s-l s + 1 s + 1/2 s-1/2 



+ ^vry + h,y^ - ^Try^ + h,y' + 0{y') 



2v^ 12 



720 



(B25) 



5 11^ B: 

^2 = - - I 2^ 



8 2\^n\in-l) J. -Ix 

^ n=2 ^ ) -) V 



dx 1 



0.440165 , 



(B26) 
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6, = 6K] 



Br, 



vn=0 

oo 



64 = 18 5^ 



n\{n — 3) 



+ 



vn=0 



(ix 1 
— 1 

" rfx 1 



n!(rz — 5) Ji — 1 x 



-9.13145- 10" 



5.9879- 10" 



APPENDIX C: CROSS SECTION 



1. Integration of the asymptotic limits 



Equation ( [4.10|) involves integrals of the following kind: 



Jo Jo 



(CI) 



where G {±z(g2(a/) ± ?2(ttj))5 ±^(?2('^/) ^ (l2i'^i))}^ J — ^'2, (see Eqs. ( [1.2|) as well as 
( pT| ) and (|D2| ) in Appendix 0) and Kj{a) = Kfj{ai,af,aci2) (see above). The asymptotic 
behavior of Eqs. (|3.19 ), ( |3.20|) , and ( |3.21|) can be summarized by the formula 



Gas{P,Zi,Z2,L,t) =g,AQ{z2- Zi)zl ''( — ) + Q{zi - Z2)zl 

I. \Zi/ \Z2 



,ifZiY-v± 



(C2) 



+ C{Qiz2-z^)zl-V- 



+ Q{z,-Z2)zl-H^-^ 



Z2 



^„|rf(zi,Z2) + a(p, Zi,Z2,L,t)Y 



The expression d{zi^ Z2) corresponds to the leading contribution C = 0. C is an abbreviation 
for the three quantities Aip~^'^^'^^ for t = and L = 00, i?i(l/^)^-'^+''ii for p = and L = 00, 
and Ci(l/L)-i+^ii for p = t = in Eqs. ( |CT| ), (|Oq ), and (g). For the quasi-Laplace 
transform of the contribution d{zi,Z2), 



V{k,,K2,L) = / dzie-^''' 

Jo Jo 



dz2e-^'''^ 



d{zi,Z2) 



(C3) 



one finds with fj = Kj/K, j = 1, 2, 



V{k^, K2, L) = P(/i, /2, LK) = 



L/l/2(/l + /2) /l/l /f/2 



(C4) 



+e 



/1/2 ^ /iVK/i + /2) 
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+ 



2l/l/2(/l + /2) 
1 



_ g-(/i+/2)LA- 



fl + /I 



/lVI(/l + /2) 



/r/2 



ln(/i//2 + 1) + (/i + f2)LK) - Et{l, hLK) 



f!f: 



J- (ln(/2//i + 1) + Etil, (A + f2)LK) - Etil, hLK) 
r / 1 V 

- 1 - Cij - In /2 - In Lir - /aLiT) 
1 ^ 



1 

hLK 



fU- 



2 

f2LK 



+ 



1 - - In A - In Lir - hLK) ]\ + 0{e 



f?f2 

In the limiting case of a semi-infinite film (L — > oo) Eq. ( |C4|) reduces to 

1 



/l/2(/l + /2) 



1 + eU 



+ ln(l + /i/A) - ln(l + ^ ) + 0{e' 



2/1 ^ 2/2 
The corresponding expression for the leading correction term 

A{h,f2,LK)= [ dzie-^^'^ [ dz2e-^'''a{p,zuZ2,L,t) 
Jo Jo 



IS 



-(f,+f2)LK 



{l + Ui + f2)LK + hULKf 



e 1 



2 flfi 

+e-/i^^(^ - 1 + (1 + hLK){CE - 1 + In /2 - InLir + /2LK)) 
+e-^2^^ ( - 1 + (1 + f2LK){CE - 1 + In /i - In Lir + /iLK)) 
+e-{/i+/2)Lxj^^^(y.^ + /2)(1 + 21nLir) 

+2(1 + InLiT) + 2hf2{LKf In Li^) | + ^(e^ 



{2 - 2Ce - In /1/2 - Ei{l, hLK) - Ei{l, hLK) 



(C5) 



(C6) 



A{h. /2, LK) = K-^+m^-^K\-^ (1 - e-^i^^(l + /iLiT) - e-^^^^(l + hLK) (C7) 



with the semi-infinite limit L — 00 



A{h, f2, LK= 00)= K-3+„^l-r,||(^/ 



^2^2 + ^^{^^ - 1 + - In /1/2 M + (^(e') > (C8) 



where C is an abbreviation for the two quantities Aip ^+''11, for t = and L = 00, and 
Bi{l / ^)^^^'^\i , for p = and L = 00. Distant wall corrections to the semi-infinite system 
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vanish exponentially. In order to obtain the analytic expressions in Eqs. ( |C^) and we 
have expanded d{zi, Z2) and a{p, Zi, Z2, L,t) in terms of e using for the e-expansion of the 
exponents r]\\ = 2 - e + 0{e^), r]± = 1 - e/2 + 0{e'^), and r] = 0{e^). The function Ei{l, z) 
is the exponential integral defined by 



00 ^-zt 



Ei{l,z)= dt = -Ei{-z). (C9) 

Ji t 

This function is numerically more suitable than the exponential integral Ei{z) 

Ei{z) = - — dt (CIO) 

J-z t 

appearing in the formulae for the correlation function. 



2. Integration of the mean-field correlation function 

Equation ( |C1|) for the full mean-field correlation function yields 

Sib = ^{p/Kf + {iK)-\ LK, /i, /2) = f\z^e-^^'- t dz^e-^'^' G{p, z,, Z2, L,t) (Cll) 

Jo Jo 

_^ 1 r 1 - e"(-^i+*)^^ 1 - e"(^i+^2)LX I _ ^-{fi+f2)LK 

= GuK '^'Yb\{f,+b){f2-b) " (f, + f2){f2-b) + (/l + /2)(/2 + &) 



e 



^{f2+b)LK _ g-(/i+/2)LX _ g-(/i+b)Li^)(l _ g-(/2+6)LX^ 



ifl-b)if2 + b) ^f^ + b)if2 + b) 

(1 _ g-(/i+fe)Li^)(l _ (,-if2-b)LK^ _ ^-(h-b)LK^^^ _ ^-(f2+b)LK^ 



e^'L^-l\ {fi + b){f2-b) + 

(1 _ e-(/i+^)^^)(l - (,-(h+b)LK^ _ ^-(f^-b)LK^^i _ ^-(f2-b)LK^ 

(/l + &)(/2 + &) (/l-&)(/2-&) 



)} 



using the notation of Subappendix |C 1| . The above formula exhibits the following limiting 
expressions: 

Sib = 0,LK = 00, /i, f2) = OuK-'^'jY-]rYrT ' (^^2) 

[Jl + J2)J1J2 
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and 



S{b = 0,LK,f,j2) = guK-'+'^{ 



1 



(/l + /2)/l/2 f?fi LK 

1 1 1 

+ 



+ 



flP2 

-_,-{h+f2)LK 



LK 



(C14) 



^(/i + /2)/i/2 flP^LK, 
Equations ( |C12| ) - ( pi4| ) lead to the following three cusp singularities: 



5(p^0,t = 0,Lir = oo,/i,/2 
5(p = 0,t = 0,Lir = oo,/i,/2: 



5(p = 0,t-^0,Lir = oo,/i,/2 



5(p = 0,t = 0,Lir = oo,/i,/2) 



= 1 



1 - 



7l + /2 
/1/2 



/1 + /2 
/1/2 



-l+r,|| 



-l+»7ll 



-l+r?|| 



(C15) 



(C16) 



and 



S{p = 0,t = 0,LK ^00, /i, h) /i/2(/i + /2) 



5(p = 0,t = 0,LK = 00, A, /2) /j2(/i + /2) 

7l/2(/l + /2) ^ 

- AVI . 



(C17) 



Lis: 



We note that the last two arguments of the nominator and denominator on the left hand 
side of Eq. ( |(J17[ ) are in general, as indicated, different from each other. For L = 00 the 
variables fj are given by — ■j(g2(a/) + q2{ai)) / K or i{q2{ai) +ql{af))/K whereas for L < 00 
the variables fj are given by —i{kq2{af) + Iq2{ai))/K or i{mq2{ai) + nq^^af))/ K with any 
combination of k, l,m,n = ±1 (see the exponentials in the last lines of Eqs. (|D1| ) and (p2|) 
in Appendix 0) . 

For the exponentiation of the p-L, ^-L, and p-C, dependences we introduce polar coordi- 
nates (see Eq. (|3.13|) ) 



CO = ^fipjKyTiLKy^, up = arctan(pL) 



i p 

= uj cos Lf, — = uj sin Lf 

LK K 



leading to the scaling function S{u!, ip, fi, /2) = S{p/K = usimp, LK 
and to its asymptotic expansion 



(C18) 



(uJCOSip) S/l,/2) 
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5(^^0,yP,/i,/2) /l/2(/l + /2 



5(^ = o,/i,/2) hhUi + h) 

_j_ g— 2 tan 



(C19) 



sin((y9)- 



-2 tan (/3 



7l/2(/l + /2) 



1 -l+'^ll 



2 ^2 
2 



Because in this section we consider only mean-field scaling functions, a simple substitution 
of the scaling variable p/K by 1/C,K in Eq. ( |C18|) leads to the same result for the ^-L 
dependences. The semi-infinite system is described by the coordinates 



UJ 



^/{^/KyTi^Kj^, <^ = arctan^ 



1 . P 

— — =CiJSm09, — =CiJCOS09 

^ K ^ 



leading to the asymptotic behavior of the scaling function 



which is independent of (p. 



fl + /2 
/1/2 



1 -'^+v\\ 



cj-i+'^il + . . . 



(C20) 



(C21) 



APPENDIX D: PRODUCTS OF WAVE FUNCTIONS 



In order to illustrate the type of transformations appearing in Eqs. ([4.10|) and ( |Cl| ) we 
present the explicit expression for the product of wave functions in Eq. ( [4.2|) : 



iq2{ai)zi 



1pf{Zi)'llJi{Zi)'llj*{z2)ij}{z2) 



(Dl) 



where s and q are defined in Eqs. ( |4.3| ) and ( [4.41 ). Thus the scattering cross section is 
proportional to a sum of 16 terms involving integrations over Zi and Z2. 
For the limiting case of a semi-infinite halfspace one has 
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V'4/2(^l)V'L/2(^l)V'S;2(^2)V'S;(^2) (D2) 

In this limit the above sum reduces to a single term. 
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FIGURES 

FIG. 1. The three scahng functions describing the lateral correlations g{pz,z/^^z/L) in the 
film (Eq. ( |3.lD ) for the limiting cases p = 0, = oo, or L = cxd: gi{u = pz) = g{pz, 0,0) 
(T = Tc,L = cx), Eq. Q), g2iv = z/0 = 5(0,^/^,0) (p = 0,L = cx), Eq. Q), and 
g3{w = z/L) = g{0,0, z/L) (p = 0,T = Tc, Eq. ( |3.4D ). The two uppermost curves correspond 
to the mean-field results for gi, i = 1,2,3, and to their feading behavior gi{xi ^ 0) = gi^i{xi), 
respectively, with xi = u, X2 = v, and X3 = w; within MET gi =52, gi,i = 92,1 = 93,1, and 53 = 53,/. 
The lower six curves correspond to 9i{xi) (Eqs. (pl|), (p^), and (p^)) and gi,i{xi) as obtained by 
perturbation flieory for d = 3. The difference between g^ and g^^i is revealed only in the inset: 
g-^ i (n) = 1 + Aiu'^+'^W , 52,i {v) = l + Biv'^+'^w , and 53,/ {w) = l + Ciu;"^+''il . Within MET one has 



Ai = Bi = Ci = -1 (Eq. (UD) and r/y = 2 whereas for {n,d) = (1,3) PT yields Ai ~ -0.9099, 
Bi ~ —1.1409, Ci ~ —0.9035, and r/y ~ 1.48. Eor vanishing scaling arguments all scaling functions 
attain 1. 

EIG. 2. G{p -^0,z,L = oo,t^ 0), G{p 0, z, L ^ cc,t = 0), and G(p = 0, z, L ^ 00, t ^ 0) 
attain their maximum value G{p = 0,z,L = 00, t = 0) = QuZ^~^ via cusplike singulari- 
ties Fi^)((p)[z(p2 + ^-2)i/2]-i+»7|i with if = arctan((pO"^), H[^\ip)[z{p^ + L-^)^/^]-^+''\\ with 
(f = arctan(pL), and Hf^\ip)[z{^^'^ + L^^)^/^]^^'^''!! with ip = arctan(L/^), respectively, inter- 
polating smoothly between the singularity Ai{pz)^^^^''\ for {t = 0, L = 00) and Bi{z / ^)^^~^^\^ 
for (p = 0,L = 00), Ci{z/L)-'^+''\\ for (p = 0, t = 0) and ^i(pz)-^+''il for (L = oo,t = 0), 
and Ci(z/L)"^+'?il for {t = 0,p = 0) and Biiz/^y^+'^W for (L = 00, p = 0), respectively. 
In 0(e) of perturbation theory (PT) the amplitude functions H[^\{p), i = 1,2,3, are given 
by Eqs. ( ^) , ( ^) , and ( ^Tl ). In ©(e) one has ^^^(0) = ijf^(7r/2) = ^1 ~ -0.9099, 
H^^\i:/2) = Hf\-K/2) = Si ~ -1.1409, and Hf\o) = Hf\o) = Ci ~ -0.9035. Within MET 
H'(\ip = 0) = H'(\ip = tt/2) = -1 and i^fV) is constant; moreover Hf\^) = Hf\ip) but not 
constant. 
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FIG. 3. The exponentiated scaling function hi{u = pz,v = z/^) (Eq. ( |3.10D ) corresponding 
to the case L = oo. We show the contour hues /ii(n, v) = h^p^i^^^.i'^ = {u^ + v'^Y^'^ , ip = arctan(t;/u)) 
for hi = 0.8,0.75,0.7,0.65,0.6,0.55,0.5,0.45 with their projections onto the uv plane as well 
as hi{u,v = 0) = gi{u) (Eq. (|3.2D) and hi{u = 0, f) = g2{v) (Eq. ( |3.3D ) which are discussed in 
Fig. 1^. The dashed lines correspond to the leading singularities gi{u ^ Q) = 1 + Aiu~^^^^\ and 
52('y — > 0) = 1 + Biv~'^^'^\'^ , respectively. 

FIG. 4. The exponentiated scaling function h2{u = pz,w = z/L) (Eq. ( |3.11D ) at bulk criti- 
cality t = 0. We show the contour lines h2{u,w) = h^p^i^ri^ — + 'uP'Y^'^-,'^ = ar:ctan{u/w)) 
for /i2 = 0.8,0.75,0.7,0.65,0.6,0.55,0.5 with their projections onto the uw plane as well as 
h2{u,w = 0) = gi{u) (Eq. {^) and /i2(n = 0,w) = 53 (it;) (Eq. (U)). The dashed lines 
correspond to the leading sing ularities c,i(n ^ 0) = 1 + Ain~^+''ll and g^iw ^ 0) = 1 + Ciw~^+'^\\ , 
respectively. In the latter case the difference between the leading behavior and the full scaling 
function g3{w) is hardly visible. Thus the leading dependence on z/L for p = 0, t = remains 
valid nearly up to the middle of the film at z/L = 0.5. 

FIG. 5. The exponentiated scaling function h3{v = z/£^,w = z/L) (Eq. ( |3.12D ) for lateral 
momentum p = 0. We show the contour lines h'j,{v, w) = h^pJiari^ — {v"^ + w^Y/'^ , 99 = arctan(w/ti;)) 
for /i3 = 0.8,0.75,0.7,0.65,0.6,0.55,0.5,0.45 with their projections onto the vw plane as well as 
h'iiv^w = 0) = g2{v) (Eq. (U)) and h^iv = 0,^) = 53(71;) (Eq. (|!|)). The dashed lines 
correspond to the leading singularities g2{v ^ = 1 + Biv^'^'^'^^^ and 53 (tf ^ 0) = 1 + Ciw^^~^^^^ , 
respectively. In the latter case the difference between the leading behavior and the full scaling 
function gsiw) is hardly visible. Thus the leading dependence on z/L for p = 0, t = remains 
valid nearly up to the middle of the film at z/L = 0.5. 
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FIG. 6. The scaling function gi{u = pz) (Eq. ( |3.2[) ) with the cuspUke singularity gi{u ^ 0) = 
1 + j4iu~"^^''II evolves out of the scaling functions hi{u,v = zj^) (Eq. ( 3.10| )) and h2{u,w = z/L) 



(Eq. ( |3.1lD ) in the limits — > and w ^ 0, respectively, which are analytic functions of u, with 
a maximum at n = 0, for v ^ or w ^ 0. The various curves correspond to vertical cuts of the 
surface shown in Fig. ^ for v = const, with w = z/L = and in Fig. ^ for w = const, with 
V = z/^ = 0, respectively. 



FIG. 7. The scaling function g2{v = z/^) (Eq. (|3.3|)) with the cusplike singularity 
g2{v — > 0) =1 + Biv~'^~^^\\ evolves out of the scaling functions hi{u = pz,v) (Eq. ( 3.10|) ) and 



h^{v,w = z/L) (Eq. ( |3.12| )) in the limits n — > and w — > 0, respectively, which are analytic 
functions of f , with a maximum at f = 0, for u ^ Q or w ^ The various curves correspond to 
vertical cuts of the surface shown in Fig. |3| for u = const, with w = z/L = and in Fig. || for 
w = const, with u = pz = 0, respectively. 



FIG. 8. The scaling function gs{w = z/L) (Eq. (3^)) with the cusplike singularity 



g^{w ^ 0) = 1 + Ciw "'^^''11 evolves out of the scaling functions h^{v = z/£,,w) (Eq. ( 3.12 )) 



and h2{u = pz,w) (Eq. ( 3.11 )) in the limits f — > and n — > 0, respectively, which are analytic 
functions of w for u ^ or w ^ 0. The various curves correspond to vertical cuts of the surface 
shown in Fig. ^ for v = const, with pz = and in Fig. ^ for u = const, with z/^ = 0. We 
note that, different than in Figs. ^ and |^, the scaling functions h^{v ^ Otw) and /i2(u 7^ 0,w) are 
nonmonotonous functions and exhibit a maximum at w ^ and a local minimum at w = 0. 
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FIG. 9. Universal scaling functions f{y) ((a), Eq. (|^)) and g{y) ((b), Eqs. (|3^ ) and 



( 3.30| )) of the film susceptibility for Dirichlet boundary conditions at both surfaces. The dashed 



lines are the mean-field (MFT) results whereas the full lines include non-Gaussian fluctuations 
obtained by perturbation theory (PT) in first order e (Eqs. ( |B22D , ( p2^ ), and p^ ). The dotted 
lines indicate the asymptotic behaviors of f{y 0), f{y — > oo), g(y — > 0), and g{y oo) given by 
Eqs. (3.40), (3.26), (3.37), and (3.32), respectively. The dotted lines correspond to the e-expansion 



of these asymptotic behaviors up to 0{e) in order to be compatible with the full scaling functions 
f{y) and g{y) whose e-expansions up to 0{e) are shown here as full lines. The dash-dotted curves 
show the exponentiated forms of the asymptotic behaviors given by Eqs. ( glo|) , ( ^^61) , (|^), and 



( 3.32| ) using the e-expansion results for the amplitudes but the best available numbers r/ = 0.031 



and ly = 0.630 for the critical exponents. f{y) has a turning point (•) at y = 1.851 in MFT and at 
y = 1.376 in 0(e); /(O) = V (Eq. (p|)). 

FIG. 10. A film (0 < z < L) filled with material 2 is sandwiched in between a 
halfspace z < filled with material 1 (typically vacuum) and a halfspace z > L filled with 
material 3 acting as a supporting substrate for the film. A plane wave with wave vector 
K* = (kj,qj) = i^*(cos Oj cos 99j, cos Qj sin sin Oj) impinges on the 1-2 interface at z = 0. The 
refiected beam has the wave vector K*" = (k*,— gj); the transmitted beam is not shown. Fluc- 
tuations in the film give rise to an off specular elastic diffuse scattering with K-^ = (kj,gj) 
= {cosaf cos If f , cos a f simpf, — sin a/), = = K. 



FIG. 11. Scaling function of the scattering cross section a (Eq. ( 4.13| )) for large pene- 



tration depths lij ^ L and vanishing lateral momentum p = as a function of the scaling 
variable y = L/^ within MFT (dashed line) and perturbation theory (full line). The dotted and 
dashed-dotted lines correspond to the asymptotic behaviors f^^^ (y) of the normalized scaling func- 
tion fo{y) = f{y)/f{0) of the total susceptibility f{y) (Eqs. (|]2|) and (|3lo|) ) in mean-field theory 



(MFT) and in perturbation theory (PT) to first order in e using in addition the best available 
exponents, respectively. 
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FIG. 12. Ratio of the normalized scaling functions of the total susceptibility (Eq. (|3.24 )) and 



of the scattering cross section (Eq. ( [4.13| )) within mean- field ((a): MET) and within perturbation 
theory ((b): PT). We use the normalization ao{y) = a{y)/a{0) and /o(y) = f{y)/f{^)- The various 
lines in (a) and (b) correspond to different penetration depths lij: I a) lij ^ L, II a) /jj ~ L, 
III a) lij S> L, and IV a) l^j imaginary (no total reflection at the interface 1-2) as marked in (b). 
The curves correspond to li = If . In the case IV a) the indicated value of L/hj corresponds to its 
imaginary part. 

EIG. 13. Scattering function S{LK; a, (3) (Eq. ( [4.18| ), full lines) and its asymptotic form 
S{LK oo;a,(3) (Eq. ( [4.22 ), dashed lines) for three different scattering geometries: Oi < Uf 



< aci2 < aci3, "i < aci3 < aci2 < a/i and < Ucu < Ocia < «/• For Ojj < aci2,ci3 the 
scattering function decreases monotonously. If one of the angles Oj or aj is larger than aci2 



oscillations emerge. This effect is enhanced if acis > Cici2- In the asymptotic form of Eq. ( 4.22 ) 
(dashed lines) there are no oscillations. In all three cases ^2,3 = 0.3- 10~^ and (a^, aj, aci2, Ocis) = 
(0.06°, 0.11°, 0.26°, 0.36°), (0.06°, 0.40°, 0.36°, 0.26°), and (0.06°, 0.40°, 0.26°, 0.36°), respectively. 
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FIG. 14. Scattering function S{p/K, LK;aij,aci2,ci3, 132,3) (Eq. (i]2|)) for t = 0. (a) shows 
the exponentiated scahng functions S{p/K, LK = oo; a, (3) (Eq. ( pl3D ) and S{p/K = 0, LK] a, f3) 
(Eq. ( C14 )) and their corresponding leading asymptotic behaviors (Eqs. ( |C15[) and ( plTj ), respec- 
tively) (dashed lines). For the leading asymptotic behavior we use the best available exponent 
r/|| ~ 1.48 and an amplitude function which is consistent with the mean- field expression (see Eq. 
( C19D ). The corrections to the leading asymptotic behavior are calculated within mean- field theory. 
For the scaling function S{p/ K, LK; a, 13) we plot contour lines {S = 0.8,0.75,0.7,0.65,0.65,0.6) 
and their projections onto the {p/K, 1/(LK)) plane (full lines) which clearly deviate from circular 
shapes, lines for ^(p/i^, Li^ = 1.5-10"'^, 3- 10-^ 4.5-10"'^, 6- 10-^ 7.5-10~'^; a, /3) (dotted lines) and 
S{p/K = 1.5•10~^3•10-^,4.5•10-^,6•10-^,7.5•10-'^,LK;a,/3) (dashed-dotted lines). In (b) and 
(c) we show the aforementioned vertical cross sections. The emergence of the {1/{LK))^^^''^\\ cusp is 
not monotonous; the vertical cross sections exhibit maxima (•) at 1/L ^ 0. The scattering param- 
eters are chosen that Ui < Uf < Ucu < Ocis with (oj, aj, aci2) Ocis) = (0.06°, 0.11°, 0.26°, 0.36°) 
and (32 = (33 = 0.3 • 10"^ 
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